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AIRCRAFT NOISE SOURCE 

AND CONTOUR ESTIMATION 

By D. G. Dunn and N. A. Peart 
Boeing Commercial Airplane Company 


1.0 SUMMARY 
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Reflecting the need for analyzing and, if possible, reducing the community noise resulting 
from aircraft operations, the Boeing Commercial Airplane Company, under contract to NASA- 
Ames has developed a computer program for predicting the noise generated by tive basic types ot 
aircraft: turbojet, turbofan. turboprop, V/STOL, and helicopter. A second program has been 
developed which calculates contours of equal noise level (footprints) and the area within the 
contours for an airplane during takeoff and approach operations. The footprint program ,s 
compatible with the NASA-Ames flight simulator. The flight simulator provides aerodynamic and 
engine performance data, and .he footprint program calculates contours for equal noise level 
thereby providing an estimate of the noise exposure produced by an aircraft operation. Typica 
results from the computer programs are shown in figure 1. These computer programs are intended 
to assist aircraft designers by identifying the noise characteristics of various aircraft and engine 
configurations. These noise levels can then be compared to community noise goals. 


Aircraft noise prediction techniques used within the aviation community are usually based on 
empirical data and the resulting procedures vary. These differences arise in numerous ways; e.g 
1 1 ) the same acoustic data can be formulated into prediction procedures with varying degrees of 
sophistication and complexity, (2) differences in noise measurements do occur in similar te= , when 
some of the important variables can not be controlled, and (3) the complexities in no.se i • oration 
and propagation have fostered more than one theoretical view of the phenomena involved. 

The source noise prediction and extrapolation techniques presented in this report represent the 
state of the art during the contract time period. The procedures are primarily empirical. Some parts 
of the procedures were obtained from published literature and in some instances unpublished 
methods used within the Boeing Commercial Airplane Company were employed. The selection o 
techniques were made to provide a base lor comparisons among aircraft design choices and lor 
evaluation of aircraft operations. However, results from these procedures can not be expected to 
agree exactly with absolute noise levels calculated by other schemes. In many instances, an 
engine/airframe configuration has its own peculiarities. These peculiarities may require corrections 


NOISE SOURCE ESTIMATION PROGRAM RESULTS 



FIGURE 1.-RESUL TS FROM NOISE SOURCE AND CONTOUR COMPUTER PROGRAMS 


Jo In- applied to the predicted levels for a particular noise source component 01 changes in one or 
more of the various calculation methods offered in this report. 

It is extremely difficult to assess the accuracy of noise prediction procedures. This difficulty 
results from insufficient data for known sources (individual component sources which combine to 
give the total noise of the system), anomalies in measured data, and the small number of different 
engine/source configurations. Comparisons of predicted with measured noise levels (Eftective 
Perceived Noise Levels. EPNdB) for cut rent airplanes have shown that the tolerance for these 
methods is generally ±5 EPNdB. It should be recognized that noise levels are logarithmic quantities, 
and an error in a noise estimate can result in a large error in a contour estimate such as area. Hence, 
the tolerance for the contour estimation procedure largely depends on the confidence level 
associated with the acoustic data that is used. 

A description of how to use the computer programs is contained in a companion report 
(ref. I ). This portion of the report contains the engineering description of the noise prediction 
procedures embodied in the computer programs. 


2.0 INTRODUCTION 

The increase in commercial aviation in the last decade has been accompanied by increased 
complaints from communities directly exposed to the higher noise levels associated with aircraft 
operations, initial attempts to reduce the community noise exposure have included changes in 
takeoff and approach procedures, development of acoustically treated inlets and mounting jet noise 
suppressors on the exhaust nozzles. Also, Federal noise regulations (ret. 2) have established noise 
limits for ih w airplanes that are significantly lower than first generation jet operation levels. This 
recent emphasis on reducing airplane noise has resulted in considerable acoustics-related research 
and development activities. These activities have been primarily directed toward defining the noise 
gem-rating mechanisms of aircraft engines and defining ways of reducing noise at its source through 
design innovations and suppression devices. 

The implementation of noise reduction technology in engine, nacelle, and acoustic lining 
design has icsultcd in a generation of quieter airplanes, e.g., the Boeing 747, Douglas DC- 10, 
Lockheed L 1011, and Cessna Citation. Numerous other programs are currently under way, each 
witli the objective of either reducing the noise of current airplanes or developing noise technology 
for application to future aircraft. Throughout these research programs, there has been only minima! 
effort devoted to developing the methodology required for predicting the total community noise 
performance of new airplanes. This report represents the state ol the art calculation procedure for 
aircraft community noise prediction. 
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The current contract. NAS2-6969. has been completed in twp parts Phases A and B (ref. 3). 
Phase A (ref. 4) consisted of providing relatively “crude” computerized procedures applying to 
noise source estimation of conventional turbojet* or turbofans. Also, part of the Phase A effort was 
the development of computerized procedures for noise contour estimation adaptable to 
flight simulation. The computer programs were designed to operate on the IBM System . / w = 

the additional requirement that the noise contour program would operate on the Xerox Sigmavu 
and VIII computers in conjunction with the NASA-Ames flight simulator. Phase B consisted of 
supplying more advanced computerized procedures for noise source prediction. The computer 
program for this purpose has been written to provide 1/3 octave band noise estimates for such 
configurations as advanced technology “quiet” engines, lift fans, lift/cruise fans, ejector/suppressor, 
blown-flap, propeller, helicopter, and tilt rotor aircraft, etc., in addition to that for conventional je 
aircraft. Also, the contour program has been updated to be applicable for the more generalized 

requirements. 


The discussion section of this report has been divided into three parts. The first deals with the 
overall view of the noise prediction procedures. The second deals with the description of the vanous 
computer modules for noise source estimation. Table 1 lists the computer modules mcluded. 


In each of the computer modules, the user has the option to specify reductions on a spectral 
basis when suppression devices are employed. For these situations where lining « •»«»“■ » 
calculation procedure applicable to optimised, single or double layer linings is avmbbte. Th«e 
procedures are included in the noise source estimation program for those items noted in the list 
the noise source computer modules in table 1 . 


After the individual noise source spectra are computed for a datum condition fftee-rield, 
1 meter from source), extrapolation technique n* then used to adjust the datum spectra to 
correspond to the noise observed as the aircraft .lies by an observer. The spectra obtained am 
functions of time and permit the calculation of the effective perceived noise level during takeoff 

and/or landing. 


The third part of the discussion deals with the procedures used to calculate noise contours. At 
the discretion of the user, an additional output of the noise source program can be a set of tabular 
data on IBM cards to define an acoustic data routine for the noise contour program. The acoustic 
function can be either peak perceived noise level or effective perceived noise level versus some 
engine performance parameter, range at the closest point of approach, and elevation angle. This 
function is then used for calculating noise contours as described in section 5.3. The computer 
programs for this task are designed to run in “real time” on the Xerox Sigma VII or VIII computers 
for flight simulations and in “batch mode” on the IBM System 360/67. 




TABLE 1. -NOISE SOURCE ESTIMATION COMPUTER MODULES 


No. of Modules 


• Measured Data 

• Jet Noise 

a. Single exhaust nozzle 

b. Co-annular exhaust nozzles 

c. Ejector/ suppressor* 

d. Slot nozzle with augmentor flap* 

e. Externally-blown flap 

• Noise Generated Inside Primary Duct 

a. Core* 

b. Turbine* 

• Compressor or Fan Noise 

a. Inlet compressor or fan* 

b. Discharge fan* 

c. Lift-fans* 

• Propeller* Helicopter and Tilt Rotor 

a. Empirical propeller procedure 

b. Theoretical propeller/rotor procedures 


Total 13 


♦Denotes optional use of lining 


3.0 CONCLUSIONS 


This is the final report for contract NAS2-6969. It describes the results for the Phase B portion 
of the contract. Empirical procedures are described which represent the srate ^t ^art an are 
best approaches readily at hand for estimating the community noise levels fo. the five bate types ot 
ulftTntioned in the summary. At the present time, comprehensive theoreucal Procures do 
not exist for noise prediction for all the aircraft mentioned. In some cases, theoretical methods can 
be computerized but they are computationally expensive, require more detailed information than is 
reaZ ava«e ."d do no, preside significantly more accurate .hrehrre ieveis than the, obtained 

by empirical means. 

Pas, experience has shown that when suppression devices are used to reduce the acoustic levels 
of major noise source(s); new noise sources always appear. New prediction P'“^ ^ 

continually developed to reflect the noise contribution of the new sources. Also the ex, sung 
procedures must be continually refined to reflect the change in engine design and the demand 

increased accuracy. 


4.0 RECOMMENDATION 


In general, aircraft noise prediction methods will change as technology improves; therefore, it 
is recommended tha, these procedures and the corresponding computer programs J 

reviewed and updated. To provide guidance on future developments and to gam a better 
understanding of the mechanisms involved, a theoretical analysis should accompany any visions 


the empirical procedures presented herein. 
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6.0 DISCUSSION 


S.l MACROSCOPIC VIEW OF TOE PROBLEM 

Noise prediction procedures have been developed for five basic types of aircraft: turbojet, 
turbofan; turboprop; V/STOL; and helicopter. Figure 2 shows these propulsion systems and the 
types of noise sources associated with each type. Of the five basic types of airplanes, the V/STOL 
aircraft is the most complicated because of the different high-hit configurations presently being 
considered-blown-flap, augmentor-wing, lift fan/jet and tilt rotor. In addition, the V/STOL 
configurations may have more than one type of powerplant, i.e., conventional turbofan wing 
mounting with lift fans in the wing or lift jets attached to the airplane’s fuselage. 

In order to provide the flexibility required in predicting the noise from the several source 
contributors for a given aircraft, the approach shown in figure 3 is used. The prediction consists of 

four steps: 

1) Solution of the Flight Path/Observer Geometry: 

In this step the airplane is assumed to move along a straight line, i.e., constant climb 
padient. At angular increments off 10°, 2© s , . . . 170° between the flight path and a line to 
the observer, sampling points for the aircraft’s position are taken whicn correspond to a 
set of observation times during the flight when the noise is heard. At each aircraft 
position, all geometrical terms required to extrapolate the noise from the source to the 
observer are determined. After the airplane/observer geometry is defined, the next step 
considers the orientation* of the noise sourceso 

2) Calculation and Summation of the Sound Levels for Each Source : 

For each point along the flight path segment, spherical coordinate angles relative to the 

noise source reference axis are calculated. See (*, * O .0 O > in flgwe4 111686 ane,6S 
describe the location on a sphere about a source where the noise is to be determined. 
Thus, non-axial-sy mme trie * * radiation patterns can be considered. 


♦The angular orientation of the poss thrust vector for a powerplant about the aircraft s lateral axis 
witti respect to the horizon. 

♦♦Except for conventional single engine aircraft, the radiation patterns will be non-axial-symmetric. 



v PRIMARY JET 

2. CORE & TURBINE 

3. COMPRESSOR (INLET) 

4. PROPELLER 


1. PRIMARY & SECONDARY JET MODIFIED FOR 
BLOWN FLAP 

2. CORE A TURBINE 

3. FAN (INLET & EXIT) 



|„ PRIMARY & SECONDARY JET MODIFIES -OR AUGMENTOR WING 

2. CORE A TURBINE 

3. FAN (INLET A EXIT) 


FIGURE 2 - REPRESENTATIVE PROPULSION SYSTEMS AND NOISE SOURCt 








1. PRIMARY A SECONDARY JET 

2. CORE A TURBINE 

3. PAN (INLET A EXIT) 

1. PRIMARY A SECONDARY JET 
t CORE A TURBINE 

3. FAN (INLET A EXIT) 

4. LIFT FAN 

JZ" y/STOL (TILT-ROTOR) 

HELICOPTER” TURBINE-OR 


Si 


1. ROTOR NOISE 

2. PRIMARY JET NOISE 

3. CORE A TURBINE 

4. COMPRESSOR (INLET) 


rmsii 


1. ROTOR NOISE 

2. PRIMARY JET NOISE 

3. COREA TURBINE 

4. COMPRESSOR (INLET) 


HELICOPTER 

INTERNAL COMBUSTION ENGINE 


ill 



LUDE 






GEOMETRIC SOLUTION of 
flight fath segment 

C«r, »»• 1 

X<* Xj, Xj, X H 

i n p w I H.L'l iJUU 

calculation and sonmation of spl*s 

FOR EACH SOURCE, £ , * X 


u engine orientation ~6g 

2) CALCULATION OF SPL*S|, 

3) SUMMATION OF SPL\ 







v 0 a VELOCITY VECTOR OF AIR SPEED RELATIVE 
TO NOZZLE 

0 * VELOCITY VECTOR FOR JET SPEED RELATIV 
J TO NOZZLE 

V * V. - V 0 » VELOCITY VECTOR OF JET SPE 
R RELATIVE TO AMBIENT AIR 


FIGURE 4.-ENGINE ORIENTA TION AND FLIGHT CONSIDER A TIONS 



The noise attributed to each source is then calculated and summed to form a total of at' 
the sound sources at each aircraft position. The noise levels thus formed represent 
free-field, index U/3 or 1/1) octave band spectra, on a sphere (radius - 1 meter) radiating 
toward the observer at each position considered along the flight path segment. 

3) Extrapolation of the Index Noise Spectra: 

In this step, sound attenuation due to spherical divergence, atmospheric absorption, an 
extra-ground-attenuation is considered (refs. 5 through 8). In addition, the interference 
phenomena of ground reflection is included as an option (section 5.1. 3.2 and 
appendix A). The results from this step represent the total (1/3 or 1/1) octave band 
spectra versus time received by the observer as the aircraft passes by. 

4) Human Response Measures: 

-he extrapolated spectra are used to calculate the human response measures; Perceived 
. Joise Level (PNL), tone-corrected PNL, and Effective Perceived Noise Level (EPNL). The 
EPNL is determined by integrating the antilogarithm (Base 10) of the tone-corrected PNL 
with respect to time (refs. 2 and 9), not from a transfer curve as was done during Phase A 
of the current contract (ref. 4). Since there is some question regarding the validity of the 
tone-correction procedure, an estimate of the effective perceived noise level EPNL is 
provided, based on the regular PNL-time history. Occasionally, the procedure gives a 
tone-correction when in fact no tones can be observed in the 1 /3-octave or narrow-band 
spectra. This estimate is denoted in the computer output by an asterisk beside the EPNL 
label. The omission of a tone penalty also solves the problem of obtaining an EPNL 
estimate when only full octave band spectra are available and tonc-corrected PNLs can 
not be calculated. Further detail on each step mentioned above is presented m the 

following sections. 


5.1.1 Definitions/Limitations/ Assumptions 
5.1 .1 .1 Flight and Weather Conditions 

The noise prediction procedures defined herein are limited to flight operations, where the 
airplane speeds are less than Mach 0.35. The SAE procedures (refs. 5 through 8) that are used for 
noise extrapolation are limited to the following weather conditions. 


Temperature 
Relative humidity 


- 1 ° to 32° C (30° to 90 F) 
30% to 100% 

0 to 1 6 km per hr ( 1 0 mph) 


Downwind 


Figure 5 shows the “weather windows” that are currently recommended by the SAE for 
acoustic testing and by the FAA for noise certification of new aircraft. 

5.1 .1 .2 Index/Free-Field Spectra 

The far-field noise data used to develop the prediction methods in this report contained 
atmospheric and ground effects that are inherent in most acoustic test data. These effects have been 
estimated and removed from the data, i.e., the data was corrected to free-fleld conditions (no 
reflecting ground plane) and extrapolated back to a distance of I meter from the source (assumed a 
point) in „'der to remove atmospheric absorption. The resulting spectra are given the term 
“Index/Free-Field Spectra." They do not represent the levels which would be observed at one meter 
from an engine, but rather far-field levels artificially synthesised in order to remove the effects 

mentioned above. 


5. 1.1. 3 Far-Field/Point Source(s) 

The acoustic far-field is defined as those distances greater than or equal to ten times the 
acoustic wavelength of interest, or ten times the characteristic source dimension. At these distances, 
the noise observed may be considered to have originated at a point. Thus, the spacing between 
engines, etc., can be considered negligible, as the observer is sufficiently far away from the airplane 
such that the noise appears to be emitting from a single point. 

5.1 .1 .4 Noise Extrapolation 

SAE procedures (ref. 5 through 8) are used for extrapolating the datum spectra (Index, 
Free-Field) to other positions in the acoustic field. The extrapolation procedures consider the 
attenuation of sound due to spherical divergence, atmospheric absorption, and the turbulent 

boundary layer near the ground. 

If the acoustic impedance of the ground is known, the interference phenomena due to ground 
reflection can be estimated (sec. 5.1.3.2 and app. At Otherwise, it is assumed that the observed 
noise levels will be typically free-field plus 3 dB -the nominal effect of ground reflection. 


5. 1.1. 5 Scaling 


For each noise component, it is assumed that the noise and the thrust from different 
powerplants can be scaled for comparisons, if the powcrplants pass equivalent mass flows when 
operating at the same gasdynamic conditions. The scale factor is determined from mass flow 


measurements: 


t 
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ACOUSTIC TESTING WEATHER: CORRECTIONS TO STANDARD DAY CONDITIONS ARE 
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STANDARD DAY CONDITIONS (REF. 6 > 

NOISE CERTIFICATION FLIGHT TESTS MAY BE CONDUCTED UNDER THESE CONDITIONS 
(REF. 2 I 


FIGURE 5.-WEATHER "WINDOWS" FOR ACOUSTIC TESTING 
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( 1 ) 


Scale Factor = / m 2 

where (ih.,iii 2 ) are the mass nows of two different powerplants operating at identical ^ namics ' 
i e Mach number and temperature. Thus, the acoustic and performance data for engine two can be 
scaled To engine one by muitiplying aU Unear dimensions, including acous^—^ 
scale factor. Note that acoustic frequency is inversely proportional to wavelength. Since thrust 
proportional to area, it scales with the square of the scale factor. 

5.1. 1.6 Multiple Engines 

For multiple engine aircraft, an increase in noise is observed over that predicted for a sugle 
engine. If N identical sources are present with no interference, the increase is 


AdB ~ 10 log 10 (N) 


(2A) 


However, it has been found (refs. 7 and 10) that the increase in noise predicted by equation <2A) is 
too high for sound which propagates near a jet exhaust of an adjacent engine in order to reach the 
observer. An empirical relationship has been developed (ref. lOl for predic.ing .he changes n 
aircraft-generated sound attributed to the attenuat, on/scattering/refract, on effects caused by jet 
effluxes. This effect is expressed as a function of the number of identical engmes N, azimuth angle 
\I/ Q , and elevation angle 0 O shown in figure 4. The formula is 


AdB 


10 


20 

N 


' COS 8 
1 + 


(1 + COT t) 


(90 c 


0.8 


) 


LOG-jq(N) 


(2B) 


This assumes that the engines have co-planar exits and that their centerlines he on a common ' 
Scan, information is available for establishing the influence of fuselage/wing shieldm . Fo 
conventional jet transports; i.e., engines mounted under the wings, the effects have 
observed. Other types of engine mountings require additional tests. 


5. 1.1. 7 Flight Effects 

The effect of aircraft motion for jet noise is accomplished by the use of the jet velocity relative 
to the ambient air (ref. 7) as the key parameter, instead of the jet velocity relative to the nozzle. 

However, an exception to this rule occurs when predicting the jet noise for an »“^“d wTth 
5 2 2 4) The overall sound pressure level data for this noise source (ref. 1 1) is normalized with 
respect to total temperature and nozzle pressure ratio. The effect of airplane speed on th, 
component is at present unknown, because part of the jet noise is generated inside the augmen o 


i 


equivalent static test, i.e., free-stream velocity equals zero. 

The other noire component (cote, turthne, fan, rotors, etc.) procedures account for the motion 
of the source by utilising the results from theory <W* »>• for the 

for flight are presented in section 5.2. 


5.1.2 Geometry Solution 

The first step in the procedure is the geometric solution for the aircraft/source p^ition versus 
time This is mquired for extrapolating the index/free-fteid spectm to the obsetver mid fo 
computin, the non-axiabsymmetric noise characteristics of each source. The analysri for 

solution is shown below: 


Required Data (see fig. 6) 

GRAD Climb gradient, i.e., tan# for Z > Z R 


X Sideline distance 

Z 0 Airplane height above the ground when at Y = 0 
Z R Observer height above the ground 

i Angle between the flight path and sound propagation path 

c A Average speed of sound over the propagation path. This value is approximated by 

whc^Czo^ndCzR arc the local speeds of sound at altitude Z 0 and Z R . 
respectively 

r. L Speed of sound at aircraft altitude, Z. This value is approximated by 

c z= ?c zo 

Mq Aircraft Mach number 
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Results (see fig. 6) 


i 


(Y, Z) Aircraft coordinates 
P Sound propagation path distance 

AP/P Relative increase in path length for ground reflected signal 

01 Angle of incidence for ground reflected signal relative to grazing incidence 

02 Elevation angle used in extra-ground-attenuation formula 

r Retarded time when sound is generated relative to the visual overhead reference; 

i.e., aircraft is at Y = 0 

t Time the observer hears the acoustic signal relative to the visual overhead 

reference 

Unit Vectors 

(i, j, k) (X, Y, Z) coordinate system 

S Direction of flight 

§ = cos0*+sin0 ic 

P^ Direction of sound propagated 

P n = -{Xi + Yj + Z n £)/P 

with Zjj = Z - Z R 

The basic approach in obtaining a solution for the flight path geometry is to solve the three 
governing equations below for the distance Y, which yields a value of Z greater than or equal to Z R . 
If the solution for Y yields a value of Z less than Z R , it is assumed that the climb gradient is zero, 
corresponding to the noise source being at the observer height. 

P cos C 3 P P N ’ S 



s ~[Y COS0 + (i 

l - Z R ) sine] 

(3A) 

P 2 

o x 2 + Y 2 + (Z 

■= z ) 2 

(3B) 

Z 

3 + Y tariQ 

0 


(3C) 
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. ' \i 


1 


A 


A value for Y is obtained by forming the quadratic equation (4), from the substitution of 
equations (31)) and <3C) into the square of equation (3A). Simplification yields 

(4) 

a y 2 - b 1 + c ■ 0 


where 


a - 

(sin|/cosd) 2 

b = 

-2 Z no tand sin 2 £ 

c = 

Z no 2 sin 2 d - (X 2 + Z no 2 ) cos 2 ! 

ii 

o 

z o* Z R 

The roots of equation (4) are then given by 

Y = C 

with 

Qi = 

b/(2a) = -0.5 Z no sin (2 d) 

q 2 = 

W-r 

Q2 = 

cos e /x 2 + (Z no cosd) 2 /tan« 


Substitution of both roots into equation (3A) eliminates the erroneous (Q, + Q 2 > rout and 
gives the solution 

u - rt _ f 

*2 (5) 


® Q 1 - Qc 


After Y is determined, Z and F are eornputed uaine equations (30 and (3B) respectively. The 
increase in path length for the ground reflected signal is computed by 

(6A) 


AP / P ■ )/ r + Y - 1 for |r| e 

for |r| < e 


O Ilw a Q 

N“* 80 ® It 

K® 1 


(6B) 
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where 


e = *7 45 
r 


D 1 

*k+l 


4 Z Z R / P‘ 
0.5 r 

he 


M» /2K - 1 \ 
■ rB <' ( ' 2K ) 


Once 
defined by 


the distances are determined. 


the angles used in the extrapolation formulas are 


0, * tan' 1 [|Z N *| /Vp 2 - (Z N *) 2 ] 


(7) 


with Z N * = Z*Z R 


P 2 - sin* 1 [ Z„ / P ] 


( 8 ) 


Next, the time that the sound is generated and observed is computed relative to the aircraft 
reference position (X, O, Z 0 ), shown in figure 6. The -esulting formulas are as follows;..e.. let 

O 9 _ ,2 


dS 2 3 Y 2 + (Z - Z Q ) 2 

dS 3 [sign of Yl l|d$ 2 
T * dS / (M Q C z ) 

t + (P / c A ) 


(9A) 


(9B) 


This completes the geometry analysis required for noise extrapolation. If the orientation angle 
(6 E ) of the noise source reference system is given, (see fig. 4), the complete geometry for no.se 

is determined using 

4 > - directivity angle 


r _fv mss. + slnfip.) / IPI ] 
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« azimuth angle 

a tan" 1 [ IX! / - (v C 0 S 6 E + Z N slnfig)] 

/t 0 a elevation angle 

. U*" 1 [ iv sln« E - Z„ C 0 SS E I / 1*1 ] 

•NOTE* 90* < *0 « '80* "hen (V «S5 E + Z„ »<«»£> 

Is positive. 

5.1.3 Noise Extrapolation 


(IOB) 

(IOC) 


The predicted noise components are treated as tree-field (index) spectra at a reference distance 
of one meter. The techniques for extrapolating from the source to the observer are basicaLy the 
same as those used during Phase A of the contract (ref. 4) with the exceptions listed below. 


1 ) Extra-ground-attenuation (refs. 7 and 8) 


2) Multiple-engine effect (sec. 5 . 1 . 1 .6) 


3) Ground reflection (refs. 14 tlurough 19) 

A revision of the extra-ground-attenuation methods was made to reflect information obtained 
from recent JT8D flight test data (unpublished) and to mateh the standards (refs. 7 an more 

closelv. 


The multiple-engine effect can cause non-axial-symmetric radiation patterns, therefore it was 
removed from the extrapolation step and included as part of the prediction methods for each noire 

component (sec. 5.2). 


A theoretical ground reflection procedure has been included as a user option (sec. 5.1 .3.2). If 
the impedance of the reflecting ground plane is known, the user of the no^e source Prediction 
computer program can estimate the effects of this interference phenomena instead of using a 
addition to the free-field spectra. 


Thus, the noise extrapolation procedures contain the following four items (see fig. 6). 


1) Spherical divergence (ref. 7) 20 logj()(P/P 0 ) 

where P 0 - 1 meter 


8 
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2) Atmospheric absorption where 3( f) 
is i he average loss coefficient 
(dB/KM) ov?r the propagation path. 
This parameter is a function ot 
frequency, ambient temperature, 
and humidity (ref. 5). 

3) Ixtra-grouiui-attenuation 

4) Ground reflection 

where (Zj/Zq) and (Ki/K(j)are 
the normalized impedance and 
wave number respectively for 
the ground. 


5(f) |P/ 1 000 1 


EGA (f, P, 0 2 ) 


GRl 


(f, P, 


AP 
P ’ 


0 


l’ Z 0 - K 0 l 


References 5 and 7 provide all necessary detail on spherical divergence and atmospheric 
absorption, respectively. Hence these items are omitted from further discussion here. 


5. 1.3.1 Extra-Ground-Attenuation (EGA) 


The best available standard, in a form useful tor prediction, is contained in reference 8. Hus 
report (based on an average of a large number of measurements) provides a procedure for 
calculating extra-ground-attenuation as a function of distance, frequency, elevation angle, and wmd 
direction. This attenuation is thought to be due to a combine non of two effects: refraction due to 
wind and temperature gradients: and dispersion due to the turbulent boundary layer. The latter is 
nearly always present; however, the former dominates for upwind propagation. Unfortunately, the 
reference defines the EGA for only a single wind velocity 10 mph. In addition, the phenomenon is 
defined onlv at source/observer heights of 1.73 meters (6 ft.) for up-wind propagation. According 
to reference 8 the attenuation is essentially constant in the downwind sector (cone angles greater 
than 120°). Because of the limitations, it has been common industry practice to use downwinc 
propagation for wind speed of 10 mph as a standard. 

Using inis standard. EGA is a function of distance, frcciucncy, and elevation angle. Data are 
shown in reference 8 at elevation angles of 0-2°. 10*. and 20° During Phase A of the contract, this 
data was represented by a function ot the following torm. 


EGAtf, I\0 2 > = 
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ft 
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where 


EGA(f, P, 0°) = 


f 


P 

*2 

K(0 


EGA at 0-2® elevation angle 

frequency 

distance 

elevation angle 

is chosen to fit the data it 10® and 20 


This function was modified during the Phase B effort because it predicted substantial amounts 
of EGA at high elevation angles, a phenomenon not observed during flight tests; e.g., nearly 2 PNdB 
at P* - 30® for 0.25 N.Mi. sideline noise estimates. During the JT8D Retrofit Feasibility Study or 
DOT/FAA, a large number of noise measurements were made at various altitudes from 122 M to 
2750 M and at several thrust values. A minimum of three flights were made at each altitude and 
thrust. From a preliminary analysis of the data, it appears that EGA approaches zero at elevation 
angles greater than 45®. In view of the above, a third formula, which goes to zero at 45° and linearly 
connects the data in reference 8 was chosen until better data becomes available. The formula is 


EGA(f. P. ft) - EGA(f. P. 0°) F(ft. f) (U) 

where EGA(f, P, 0®) is obtained by linear interpolation with respect to log(P) and log(f) on the data 
given in table 2. The function F(/K 0 is shown in [ ' &lK 7 antl ,s tabu,ated ,n tab,e 

5. 1.3.2 Ground Reflection 

The particular model considered is a point source, homogeneous media (air and ground), and a 
smooth/inflnite/reflecting plane with complex acoustic wave impedance based on the acoustic 
analog to that in electromagnetic theory (refs. 14 through 19). From tins model, it has been tound 
that the reflection effects are quite sensitive to Ihe aircraft/observer geometry, the source 
frequency, and wave number ratio (K,/K 0 ). end the normal impedanee (Z,/Z 0 ) ol the ground. 
Complication aria’s because the values of these last two parameters are not common knowledge lor 
various types of terrain. If the parameters Z,/Z n and K,/K 0 am not known but only guessed, he 
predicted ground reflection effect can make the Effective Perceived Noise Level more in error than 
that obtained by simply adding 3 dB to free-ficld data. However, this phenomena does deserve 
study because its spectral effects are significant. A detailed analysis bar the optional ground 
reflection procedure incorporated in the computer program is provided in appendix A. 
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FIGURE 7 -DECAY FACTORS FOR COMPUTING EXTRA GROUND ATTENUATION 
(COMMERCIAL OCTAVES ) 




TABLE 3.-DATA POINTS FOR EGA DECAY FACTOR 



COMMERCIAL OCTAVE BAND NUMBER 
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Values in Table above were obtained from Figure 
4 (Appendix 2) of SAE AIR 923. The AdB error 
was minimized by using 

f = 2J * AdB J w / £ Acl& j! 

for ranges P = 250, 350. 500, 700, 1000, l-»00, 
2000, 2800, 4000 feet. 










5.1.4 Lining Treatment 


Current technology in predicting the attenuation of acoustic linings incorporates a combina- 
lion of experimental correlation and theoretical analysis. The acoustic wave attenuation analysis is 
based on a rectangular duct with mean flow and boundary layer effects. Equivalent duct heights for 
non-rectangular duct geometries are obtained by equating flow areas, treated areas, and duct 
lengths. Far-fleld attenuation directivity corrections have been obtained from engine ground 

test data. 

In the application of this technology, three types of prediction procedures of different detail 
and complexity can be identified. These three types are (1) prediction of the attenuation spectrum 
for a given lining design, ( 2 ) prediction of the lining parameters for a given duct configuration where 
the attenuation of a noise spectrum is maximized, and (3) prediction of the attenuation spectrum 
for a given duct with lining parameters unspecified, but assumed to be chosen such that t e 
attenuation of a given noise spectrum is maximized. 

The first kind of prediction (l)can be made relatively inexpensive by the solution of the 
equations governing wave propagation in a lined duct of somewhat idealized geometry leading to 
solutions which compare reasonably well with data. To accomplish the second kind of prediction 
requires an optimization program. The optimization program iterates the procedures contained in 
the first kind of prediction (1) resulting in the optimum lining parameters that maximize the 
attenuation of a noise spectrum. Although optimization programs exist, they are too costly to run, 
except in final design of a lined duct configuration. 

The greatest need for predictions are of the third kind, which arise in trade studies, where the 
effects of such parameters as inlet length, number of splitter rings, and engine choice are 
investigated. For these cases, a simplified procedure which uses attenuation spectra corresponding 
to optimized linings is used. These attenuation spectra are somewhat idealized, i.e., made to 
confonn to a standard shape in order to avoid the expense of an exact calculation. Experience with 
this approach, though approximate, shows good correlation with results from detailed predictions 
for perceived noise level reductions. This approach has been incorporated into the computer 
program. The following discussion pertains to the type 3 prediction procedure. 

This procedure encompasses two types of optimized linings: single-layer and double-layer. The 
single-layer procedure gives the user two options. The first option considers only a single design 
point; i.e., fixed engine conditions. The second option considers multiple-design-points, i.e., the 
engine conditions vary over a limited range. A compromise between peak attenuation and 
bandwidth is made because the lining is expected to attenuate tones that track with the engines 

shaft speed (rpm). 


For double-layer linings, an increase in bandwidth can be realized with the same peak 
attenuation as that given by single-layer linings. Thus, double-layer linings can be used for the 
multiple-design-point option described above, to provide the same perceived noise level reduction 

for a slightly less amount of lining «.ea. 

The source noise computer program contains lining attenuation calculation procedures as an 
option for lining treatment of the following noise components : 

1 ) Compressor and inlet fan 

2) Discharge fan 

3) lift fan 

4) Core and turbine 

5) Ejector-suppressor jet noise 

6) Slot nozzle with augmentor-flap jet noise 

Within these procedures, there are several methods available for calculating the lining 
attenuation spectra. These methods are as follows. 

1) For each target frequency, the user defines the magnitude of maximum attenuation and 
the percentage of the total area that is treated. The program then determines the 

spectrum shape. 

2) The user defines the effective duct height, the ratio of apparent treatment length to 
effective duct height, and the percentage of the total area that is treated for each target 
frequency. The program then determines the spectrum shape. 

3) The user defines the geometry of the lining in terms of the length and radii of cylindrical 
walls, and the percentage of the total area that is treated for each target frequency. The 
program then determines the spectrum shape. 

The user is limited to the configurations shown in figure 8 when defining the linings 
geometrically. 
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a) CIRCULAR DUCT 

(Number of walls - 1) 
(H * R) 


b) ANNULAR DUCT 

(Number of walls * 2) 
Innermost and outermost 
walls are lined on one 
side. 

(H = R-j “ ^ 


c) "n" CONCENTRIC WALLS 

(Number of walls a n) 
Innermost and outermost 
walls are lined on one 


The lining attenuation prediction precede involves the successive use of ftve figu'ea. Th^e 
curves consider duct geometry, treatment area, target frequency, spee o soun , 
number, attenuation ow* a ran* of power settings, attenuation specn.m : 
angle. Depending on user requirements, any or all of the procedure may be used. The use 

figure is explained in detail below. 

Figure 9 is an estimate of the peak attenuation obtainable for an optimum singWayer lining at 
a single power setting and a zero Mach number. The required data are: 

n L/H One half the ratio of actual treatment area to duct cross sectional flow area. 

* * ' . . |, K a fifzoL n f that which would be calculated 

Actual treatment area is typically about 65% ot tnat wmcn 

from a nacelle halfsection drawing. L is the apparent treatment length. H is the 
effective duct height. See sketch for example. 


Inside of outer cylinder Is lined 
Outside of inner cylinder Is lined 


» f t H/c Non-dimensional target frequency, where f t is the target frequency for peak 
attenuation, and c is the speed of sound. 

Figure 10 shows the variation of the single-design-point peak attenuation with ‘he ductMach 
number However it should be remembered that these curves represent optimum linings at the 
ZTv'cy and different Mach numbers; no, the same lining a, a different Mach number, as rs the 

case in typical duct data and theoretical am iysis. 

Figure 1 1 shows the compromise when a lining is designed to operate effectively over a range 
811 • *u i nop The correction factors cause a reduction 1C in peak attenuation, 

of power settings-the usual case. The respectively for single or double layer 

UnL^ m0re than ^ eXhaUSl m<>de beCaU t e ’ 

ring's changed, the engine blade passage frequences end the ££*£££ 
frequency shift in opporite directions tar the inlet mode, and m .he ssme drmetron for 

exhaust mode. 

Figure 12 gives the attenuation spectrum shape. This is used to obtain the attenuatim. of 
sound a, frequencies other than the target frequency. It is considered representative of an average 
single-layer-lining ease: the reasonable corpromise between the maximum obtamable attenuation 
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FIGURE 9.-EFFECT OF DUCT GEOMETRY, TREATED AREA, TARGET FREQUENCY, AND 
SPEED OF SOUND ON REALIZABLE LINING ATTENUATION 


k dB) M /(AdB) 0 FOR SINGLE DESIGN POINT 



II- MACH NUMBER 

( AdBfc - REALIZABLE ATTENU ATI ON 0 M 
(AdBf o " REALIZABLE ATTENUATION 9 M - 0 
(AdB)n-(AdB>o [<AdB ^/<AdB > 0 ] XK 

K - MULTIPLE DESIGN POINT CORRECTION ( FIGURE 11 ) 





MACH NUMBER 


INLET'' 


EXHAUST 


FIGURE 10.— ATTENUATION CORRECTION FOR MACH NUMBER EFFECT 





SINGLE DESIGN 
POINT OPTION 


MULTIPLE DESIGN 
^ POINT OPTION 


(M ) 


-OJ 0 0.5 1.0 


OPTIMIZED- SINGLE- LAYER LININGS 


INLET MODE 


EXHAUST MODE 


( KM 



NOTE 

USE K - 1.0 FOR AdB 
CORRECTION FORMULA 
IN FIGURE 10 


(M ) 


-0.5 0 0.5 1.0 


OPTIMIZED- DOUBLE-LAYER LININGS 


FIGURE 11. -BANDWIDTH CORRECTION FACTORS 




FIGURE 12.— TYPICAL SHAPE OF LINING ATTENUATION SPECTRA 
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the target frequency and the requtrement for attenuation bandwidth using the stngk^gn-potn. 
ZL TO appr.tdn.ate the use of the utuitipMesign-poin. option or the use of doubie .aye 
linings, the bandwidth is increased by multiplying the argument in octaves. log 2 (f/f T ), by the factor 
(K or K'), shown in figure 1 1 . 

Figure 13 gives the far field directivity correction to apply to the attenuation spectrum. Figure 
14 shows a comparison between the predicted attenuation and the experimental data. 

5.1.5 Configuration Corrections 

The effects of fuselage/wing shielding and reflection have not been defined due to the scarcity 
of information on the influence-of engine placement and, furthermore, it is expected that these 
effects are small for conventional jet transports. However, any radical change in engine loca .on, 
e.g„ oventhe-wing-mounting, or use of suppression devices, etc., could result in a subst.nt.al ctange 
in noise level and require corrections to be employed. Since .here is no way to 
changes in airplane/engine configurations, the approach taken here is to let the user of the prognun 

define corrections for each noise component. 

The corrections prescribed by the user for each noise component can be used for: (1 > a A dB 
to be subtracted from the overall sound pressure level or (2) AdB’s to be subtracted from the "owe 
spectrum. Hie prognun permits the corrections to vary with .he direclivity angle, ♦ . he 
Doppler-shif. option is selected in predicting a noise component, .he program assumes he 
corcections are representative of that obtained from a static test and it will Dopplenshift 
prescribed correction spectra. 

5.1.6 Summation of Component Noise 

T1.C problem is stroctured to permit the calculation of the datum spectra for each noire 

component on a common reference sphere at angles HP. 20". 170*relati.e to .he fltgh. path. 

Thus, individual noise components are added together in the usual manner for logarithnuc 
quantities. As each noise component is determined, the total noise for all the sources ri accumulated 
using equation (12) below. The result represents the total sound radiating toward the observer at 
each aircraft position considered. Further detail on the prediction of the component no.se fo. each 

source is given in section 5.2. 
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3.1.7 Output for Noise Contour Estimation 

One of the requirements of the noise prediction computer program is the iinkage of its output, 
a data table, with the noise contour program. The output variables for the data table are listed 
below and illustrated in figure 1 5. The normal used for the output file is (1PE12.3, 3E1 2.3). 

ML Noise level (EPNL or PNL maximum) corresponding to each district combination of 
values for EPP,® , LR. 

EPP Engine performance parameter. This variable is to be specified by the user as a correlation 
parameter for noise and may correspond to engine pressure ratio, jet velocity, engine 
speed (rpm), etc. During Phase A of the contract, this parameter was engine pressure ratio 

(ref. 4). 

a Elevation angle in degrees (ref. 4). It is computed as l o = cos' 1 (X/P) 1 where P is the range 
at the closest point of approach (fig. 6). 

LR Logarithm (base 10) of the range at the closest point of approach (CPA). It is computed 
as logjo(P) when f is equal to 90° (fig. 6). 


5.2 NOISE SOURCE ESTIMATION 

The computer program has been generalized to accomodate several difterent types of noise 
sources associated with current and future aircraft. The following list describes the types of noise 
source prediction modules that have been included in the “Noise Source Computer Program. 

Reference 


1 ) Measured Data 


2) Jet Noise 

a) Single exhaust nozzle 

b) Co-annular exhaust nozzles 

c) Ejector/Suppressor nozzle (JT8D design) 

d) Slot nozzle with augmentor flap (STOL) 

e) Externally-blown flap (STOL) 

3) Core and Turbine Noise (Turbojets/Turbofans) 


7, 20. and 21 

22 

23, 24, and 25 
1 1 . 26, and 27 
28 and 29 


38 


VARIABLES 


) ML 

> EPP 
I) LR 
DO 

♦MOTE# 


A THREE-DIMENSIONAL DATA ARRAY OF NOISE LEVEL 
AS A FUNCTION OF(EPP, LR.O t THE NOISE LEVEL 
VALUES REPRESENT EPNL OR PEAK PNL, ETC. 

A ONE-DIMENSIONAL DATA ARRAY OF ENGINE PERFORMANCE 
PARAMETER VALUES FOR THE <NL> ARRAY. 

A ONE-DIMENSIONAL DATA ARRAY OF LQG^ RANGE AT CPA) 
VALUES FOR THE (ML) ARRAY. 

A ONE-DIMENSIONAL DATA ARRAY OF ELEVATION ANGLES 


FOR THE (NL) ARRAY. 

,»T* ARRAYS ML. EPP. LR. O 1 DEFINE THE THREE-NNNENSIONAL 
IABULAR FUNCTION. NL ■ 'x<EPP. LR. « > »« CORRESPONDS 

0 LEVEL FLIGHT AT A NOMINAL AIRCRAFT VELOCITY, AND 
MormuiTv Aufti F OF PEAK NOISE RADIATION, 




4) Compressor and Fan Noise 

a) Conventional turbojets/ turbofans 

b) Lift-fans 

5) Propeller, Helicopter ana Tilt Rotor Noise 

a) Empirical propeller procedures 

b) Theoretical propeller/rotor procedures 

All, or any subset, of the modules available can be used for noise prediction. Each module may 
be called up to a maximum of three times, corresponding to an aircraft having three different types 
of engines mounted on it. This restriction does not apply to the number of engines of the same type 
and orientation angle, 6jj. 

In each noise prediction module, the known non-axial-symmetric characteristics in the 
radiation patterns are considered in the calculation of the datum spectra (free-field, index); e.g., the 
multiple-engine correction formula equation (28) in section 5 .1.1.6. This also includes any user 
specified installation effects (sec. 5.1.4 and 5.1.5) on the radiation patterns. Other factors 
considered are flight versus static conditions; relative velocity, Doppler-shift, etc., (sec. 5. 1.1. 7 and 
fig. 4). The optional use of lining as a suppression device can be included in the modules; 2c, 2d, 3, 

and 4 listed above. 

As evident from the options available, the user of the program can predict the noise for almost 
any aircraft presently flying, and some which have not been built yet. However, empirical 
procedures have their limits and failure is anticipated. (The complexity of the types of aircraft 
presently being studied could also be the limit.) For these cases, a module has been included to 
accept measured data to provide more accurate results when one or more of the models presented 
are considered inadequate. 


30 and 31 

32 

33 through 38 


5.2.1 Measured Noise Data 

Since measured data is generally more reliable than that obtained from current prediction 
procedures, the capability of including measured data in the computer program is available. The 
measured noise is assumed to be sound pressure level spectra, SPLS, given in dB re 20PN/M as a 
function of three or four variables for axial-symmetric or non-axial-symmetric type sources, 
respectively. The independent variables are frequency (eight preferred 1/1 octave bands or 
twenty-four 1/3 octave bands defined in table 4), some prescribed engine performance parameter, 
directivity angle (*). and elevation angle ( 0 O ). Prescribed spectra are assumed to represent far-field 
noise extrapolated back to a free-field, index condition. Thus, the levels can be treated is 
independent of local ambient conditions. 
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TABLE 4 -FIL TER BAND DEFINITION AND ATMOSPHERIC ABSORPTION 
AT 45.7 m ( 150 FT I FROM A SOUND SOURCE 


TYPE OF 
ANALYSIS 


1. PREFERRED 
FULL 
OCTAVES 


17-19 

20-22 

23-25 

26-28 

29-31 

32-34 

35-37 

38-40 


CENTERS 

FREQUENCY 

f 

(Hz) 

63.1 

126. 

251. 

501. 

1000 

2000 

3980 

7940 


FILTER 

LIMITS 


44.7/89.1 

89.1/178 

178/355 

355/708 

708/1410 

1410/2820 

2820/5620 

5620/11200 


ABSORPTION \ 
COEFFICIENT 

ofc(f) 

(dB/KM) 


ATMOSPHERIC! 

absorption 

Am 

(dB) 


preferred 

(1/3) 

OCTAVES 


17 

50.1 

18 

63.1 

19 

79.4 

20 

100. 

21 

126. 

22 

158. 

23 

200. 

24 

251. 

25 

316. 

26 

398. 

27 

501. 

28 

631. 

29 

794. 

30 

1000 

31 

1260 

32 

1580 

33 

2000 

34 

2510 

35 

3160 

36 

3980 

37 

5010 

38 

6310 

39 

7940 

40 

10000 


44.7/56.2 

56.2/70.8 

7C.8/89.1 

89.1/112 

112/141 

141/178 

178/224 

224/282 

282/355 

355/447 

447/562 

562/708 

708/891 

891/1120 

1120/1410 

1410/1780 

1780/2240 

2240/2820 

2820/3550 

3550/4470 

4470/5620 

5620/7080 

7080/0910 

8910/11200 


[i> 15°C, 70% Relative Humidity 
h > 3(f) x 0.0457 KM 

&> 
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The computer program’s function is to interpolate on these data at specified aircraft operating 
conditions and extrapolate the spectra to the observer. No effort is made to scale, correct or 
aircraft speed, etc., because no information is known about the type of sound source the noise 
represents or the measurement conditions in which the data were obtained. 

5.2.2 Jet Noise 


For the procedures presented in this report, jet noise is defined ss the noise generated by jet 
flows as they exhaust into the atmosphere. The actual noise generation is thought to take place tn 
the flow regions where the jet flow interacts with the atmosphere. The noise generated upstream 
inside the engine is discussed in sections 5.2.3 through 5.2 A4 , 

In the case of a single jet, the noise producing regions are shown shae in figure 16. In the 
past the noise produced by this jet was correlated with three parameters-denstty, area, and 
velodty relative to the ambient air, yielding considerable success in predicting maximum passby 
noise of turbojet engines (refs. 4 and 7). 


However, turbofan engines have replaced the turbojet engines as tne most common power 
plant for todays air transport fleet. The new turbofan engines are considerably quieter than e 
turbojets and the jet noise produced by the newer engines has been observed to differ from tha 
predicted by the SAE procedure (ref. 7). Clearly, a revised procedure was necessary. Dunng P ase 
of the current contract, a relatively “elude" revision was made (ref. 4), but the subject demanded 
better precision for estimating time-integrated subjective measures, as prescnbed for the Phase B 
part of the contract. This required prediction of noise at several angular positions relative to he 
inlet centerline, instead of just at that angle where the maximum noise occurs. For the case of lift 
jets attached to an aircrafts fuselage, it was of particular importance to calculate the relative jet 
velocity vectorally as shown in figure 4, section 5.1 because an increase in jet noise occurs wi 
crossflow imposed on the jet, i.e., the relative jet velocity is greater than that without crossflow. 
The procedures presented here for jet noise are empirical and represent the state-of-art m solving 

this problem. 


Before the noise prediction procedure for a single jet Is presented, it is worth noting some 
technological developments (rer. 39 through 41) that are taking place at this time which could 
result in a more comprehensive prediction procedure. Figure 16 shows the noise producing shear 
regions in a single jet and the corresponding relative acoustic power level as a function of distance 
from the nozzle exit plane. This distribution of energy can be further broken down with respec o 
frequency f/f„ and position X/X 0 as was done in references 39 and 40. Hence the noise P™*“«« 
regions appear to be structured and not random as previously believed. Thus, it may be possible that 
a one-dimensional source distribution model can be developed for predicting the no.se in a far-field. 


42 




This math model could then be used to cornet ground static noise data to a point source, and 
,he parallax errors that occur in extrapolating the data to oth,. positions in the accuse 
field The model could also be used to simulate flight conditions from ground static tes'i by 
accounting for the change in potential core, velocity, etc., that cause a spectrum shift with flight 
speed. The present SAE practice (ref. 7) gives two spectrum shape curves f,-r flight and ground 
static predictions when in theory there should only be one. 

In addition to the above, a better normalization of jet noise results when the acoustic power is 
related to the mechanical power, convection Mach number, density and temperature ratios as 
outlined below in equation form and illustrated in figure 17. 

Total acoustic power. 



where the one-dimensional jet flow parameters are 


Wo = 

mechanical power - (^) 1 Vj - Vq| ^ 

It 

mean density ratio 

sV 

II 

mean static temperature ratio 

M r = 

mean convection Mach number 

Stf 

0.5 |V,-V 0 |/C, 

m = 

5 for quadrupole sources 

A 

Vj - 

jet velocity vector relative to the nozzle 

a* 

1 ! 

nozzle velocity vector relative to ambient air 

Cj ■ 

mean speed of sound in the jet 

m - 

mass flow = PA J Vj 1 


m 


«pw 1 c* A Vl 01901 01 “2i2£ 


DATA FROM FIGURE SC OF REFEREMCE 20 


SLOPE 


n- h 3 SLOPE 

(UPSTREAM FACILITY NOISE) 


REFERENCE CONDITIONS: 

-j- DENSITY - 16.02 KG/M 3 ” 1 LBM/ft 3 
AREA- 0.0929 M 2 -lft 2 
TBIP.(T S0 ) • 298.2* K - 534.7* R 
DISTANCE (R> - 0J04B M - 1 ft. 

VELOCITY - 0.3048 mps - 1 Ip* “ 

STATIC TEST: - 0 


MACH NO. ~ V|/Cj 

FIGURE 1 /.-SPACE-AVERAGED OVERALL SOUND PRESSURE LEVEL ASA 
FUNCTION OF MACH NUMBER 











Equivalent forms of equation (13) are: 


W, 


~ W ° (£o ) ( Tso) 


I. 5 


V j-Vq 

Co 


OR 




-v. 


The above suggestions are preliminary and tether work is required prior to r.. rporation into 
the empirical procedures currently used. This should not be interpreted to imply that the methods 
now used are inadequate and cause gross emus, but it does suggest that they will reqinre additional 

refinement to increase the range of application. 

5.2.2. 1 Single Exhaust Nozzle 

In the present procedure (ref. 20), the noise produced by the jet of a single nozzle has been 
correlated with relative jet velocity (0, - V 0 ) as shown in figure 18 . This curve was developed from 
the data given in reference 20 for a 2.54 cm diameter nozzle and was extended to match the 
predicted results of reference 7 for velocities greater than 760 mps (2500 fps). In formula form, uie 
space-average, overall sound pressure level is 


mfi ~ F(V r ) + 10 Log 10 (tt-) 0 


(14) 


where 

P is the jet density and the reference density, Pr, is 16.02 KG/M 3 (1 Ibm/ft ) 

A is the fully-expanded discharge area and the reference area, A .» ' 0.0929 M 2 (1 ft ) 

V R is the magnitude of (Vj - V 0 ) 
n varies with V R as sho wn in figure 19 

Previous definitions (ref. 7) for the overall SPL used a P 2 normalization instead of aP" term 
used here. The value of n was determined by “force-fitting” the formula to expenmantal data (ref. 
20) so that the term, F(V R ), was approximately proportional to V R 8 for velocities less than 

760 mps. 
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The free-fifild, space-average, SPL spectrum for (1/3) octaves is determined from figure 20 as a 
function of Strouhal Number. 


S?L(f) = OTFT + F(SN) 


(15) 


' where 

f ■ geometric-mean frequency for a pass band 

| SN = f/f Q 

f 0 * characteristic frequency = Vr/D q 

Do I 2 - 

I . 

Two spectrum shaping curves have been provided in figure 20 corresponding to the flight and 

ground curves given by SAE (ref. 7). 

Next, the SPL spectrum at a particular directivity angle ( 4 >) is obtained by interpolating from 
* the data shown in table 5. This table was developed from data given in reference 21 and provides a 

correction to add to the space-average SPL for the desired spectrum. 

SPU(f) - sPT(f) ♦ F(S. v„, I/O < 16) 

I 

where 

S = an effective Strouhal number to enter the table 

t - ftv u o 

u 0 - 304.8 m/s » 1000 fps 

The results from equation (16) represent the free-field levels for a sngle engine at 45.7 M 
I (150 ft) from the source. The spectrum is corrected to the datum condition (P = 1 M) through use 

of table 4. . i 

SPL(f) = SPL(f) | + 33.2 +AdB(f) 

\ M 45.7 M (17) 

I The effect of multiple engines of the same type and orientation is accomplished by adding the 

correction (eq. (2B), sec. 5.1. 1.6) to the result given by equation (17). This completes the 
prediction of jet noise for a single nozzle. 
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S.2.2.2 Co-Annular Exhaust Nozzles 
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For a co-anmilar jet, i.e., a jet with both primary and secondary flows, the norse produci g 
legions are shown In figure 21 and am defined as: (l)the inner shear layer, which rs due to he 
interaction between the prime* and secondary flows; (2) the outer shear layer, which .s dueto he 
interaction of the second.* with the surrounding air; and (3) the runted flow regron where the 
combined jet flows have become fully developed turbulence. 

For noise prediction purposes, the ctwnnular jet is considered to have the noise 
characteristics of two independent single jet flows; one represented by the rimer shear layer and 
second represented by the summation of the outer shear layer and the mixed flow region. 

The noise characteristics of the inner shear layer alone are predicted first as though _<he 
seconds* flow was absent, i.e., the srnne as in section 5.2.2.1, except that the levels am then 
adjusted to account for the presence of the seconds* flow surroundmg the primary flow. T 
predict the noise of the outer shear layer and mixed flow region, it ts necessa* to calculate tte 
acoustical equivalent flow parameters for the second.* stream of the co-annular flow sys *' m ' 
is a phenomenological, force-flt approach and it has no physical rmphcatrons to the m ean 
one-dimensional flow parameters for the flow region being considered. These c«^pmmnet^ 
am then used to predict the second.* flow jet noise as though it was a srngle jet. The predicted 
Tise of the total ^annular jet is the energy sum o, that produmd by the two flows mentroned 

above (ref. 22). 

Noise prediction for inner shear layer. -The noise from the inner shear layer is predicted in the 
same manner as that described in the pmvious section 5.2.2.1 with the folkrwrng excep ro . In the 
step where the space-averaged SPl spectrum is calculated, (eq. 15), a correction term rs mserted 


ffl = SPL(f) 


AdB(f) 


(18) 


where 


AdB(f) = 


m 


0l’ V J2 " 


10 logjo 


VJ1- V J2 

Vjl 


m 


experimentally determined exponent (ref. 22) shown in figure 22 
F(A2/Aj, f/f | ) 

primary/secondary velocities relative to the nozzle 


I 






Aj ( A '2 = primary/secondary discharge areas 

f, = Vjj/Dj 

D, 2 - T A 1 

Note: A singularity exists whenever V J2 approaches Vjj. In the program whenever |(Vjj - 
Vm )/V,i I is less than 0.1, the inner shear layer is assumed to vanish and hence produces no noise. 
Further, the above procedure has not beer verified for the case of V J2 > Vj|, but the computer 
program still considers this a valid case. 

Noise prediction for outer shear layer and mixed flow region. -In this step, the noise is 
predicted in the same manner as that described for a single exhaust, section 5.22.1, except the 
parameters (P, A, V R ) are the acoustical equivalent terms defined below. 

p = mean one-dimensional flow density of secondary discharge 


A = A, + A 2 



Vr 2= Vj 2 - 2 V, V Q cos «+ V Q 2 


a = angle betwe en gross thrust vector and the direction of motion 

5.2.2.3 Ejector/Suppressor Nozzles 

A multi-element suppressor nozzle is shown in figure 23. This modification (lobe or tubular 
nozzles) of the exit hardware of jet engines can yield a considerable amount of noise suppression 
when compared to a conventional circular discharge nozzle. The suppression is believed to result 
from the change in turbulent mixing-an alternation of the turbulence scale and a reduction m the 
mean relative jet velocity gradients (ref. 24), since an increase in induced secondary air is observed. 
This implies an integration over the total volume of the jet. However, this is just a gross observation 
of a very complicated phenomena. In fact, the processes are so complicated that they defy 
theoretical analysis. Only empirical methods have yielded feasible designs (ref. 25). 





For the purpose of jet noise prediction, the noise from e multi-element nozzle can be 
--.-w. to consist of two parts-(l) premerging noise and (2) postmerging noise The jn.wnerg.ng 
noised generated the region dose to the nozzie where the structure of the indmdual je s 
identified. The postmerging noise is generated in a region downstream ta*. £ 
individual jets and induced secondary air has merged into a sutgle um 

velocity The high-frequency portion of the resultant total jet noise spectrum is usually dominated 
by the ptemerging noise, while that of the low-frequency portion is assorted with the postmerpng 
noise. The noise for each component is predicted in a manner similar to that of a round nozzle^ The 
total jet noise is then obtained by summing, on an energy bash, the spectra for premergmg 

postmer^jng noise. 

When a shroud, commonly called an ejector, is added to the multi-element suppr^. an 

, in suppression can occur, provided that the shroud length to diameter ratio, L/D, is larg . 

However long ejectors have considerable weight, and losses in flight associated with them; hence, 
“e not hl studied in depth to a nohe suppression item. Another approach (ref. 25) toes a 
shorter ejector which incorporates lining to achieve the same result. Even t s a PP roac * * 
because all the premerging noise does not propagate normal to the ejector walls, and 0 y p 
rrhel in.eroep.ed by 8 , he lining. Furtiter, only limited types of lining materials can be used due 

to the thermal environment of the exhaust. 

For short shrouds, L/D less than 2.S, without lining, no significant reduction in premergins 
noise! been obsereed when compared to the, of a' “bare” suppressor configuration. Thh ieato to 

the pHon tha, the noise from a short, hardwall ejector/suppressor can be predated m 

manner similar to that for a “bare” suppressor. This requires knowledge of the ejector P erfor ™ n “’ 
however, because the presence of a shroud (ejector) imposes a constraint on the boundaries o 
premerging and postmerging regions. Ejector performance can be obtained by use of a th eoretical, 
one-dimensional, flow analysis if one assumes 100% mixing inside ihe shroud. Appendix Bco 
an example of the parametric curves that can be obtained ftom the configuration shown m figure 

24, taking this approach. 

The present noise prediction procedure is essentially empirical. It is based on the analysis of an 
extensive amount of round nozzle and suppressor noise date from the fofiowmg types of tests. 

1) Full scale JT8D, JT3C, JT4/J75 and JT12 static engine tests. 

2) Model scale hot flow test (A -45 .6 cm ). 


3) Flight test for the 707, 727, and 737 airplanes. 


ENGINE AIR INLET 




I 


The results of the procedure are 1/3 octave band spectra at the frce-field, ,nd ^ 
described in section 5.1.12. The overall accuracy of the procedure is dependent upon the ability to 
predict the premerging noise spectrum. This, in turn, depends upon the 

"pZX upstream from the suppressor nor* exit which 

increasesThe turbulence in the jet for the JT8D configuration and hence increases the P' cmcr ^ 
“he lts. in table 6 represent the tolerances in PNdB based on observa.rons from ^ndriahc 
tests that can be expected for the procedure, provided that the area ratio* lor the «PP»«« ' » 
between 1 ,5 and 4.5. The area ratio is defined as total flow area (primary plus mduced seconds 
divided by the primary discharge area at the suppressor exit plane. 

Figure 25 shows a comparison of a sample prediction with measured data. 

Postmerging noise prediction. -Consider the ejector/suppressor configuration shown in figure 
24 The postmerging noise for the ejector exhaust is assumed to be simdar to that o. a conven ton 
M techmques described in secton 5.2.2.1 could be applied; however a shghfiy 
different approach is taken here. The overall sound pressure level for a single engmeri reWed toj the 
relative jet velocity, density of the exhaust, stafic temperature and dtscharge area. The relatton 

OASPLWO * f,(v r . f ) + io (-sj <l9) 


where 


F,(V R , 4> ) is obtained from figure 26. 


relative jet velocity = jVj - V Q 


A A 


Vi = mean one 


-dimensional flow velocity for the ejector exhaust 


p = mean one 


dimensional flow density for the ejector exhauct 


reference density - 16.02 KG/M^ (1 lbm/ft *) 


Te = mean one-dimensional flow static temperature 
^ SR = reference temperature from figure 27 
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TABLE 6-TOLERANCES 



a Short. hardwall shrouds: L/D 2.5. 
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FREQUENCY IN HERTZ 

NUMBER OF TUBES 37 TOTAL TEMPERATURE 703°l( (1205 8 R) 

AREA RATIO 3.3 AREA 0.352 M* (3.78 FT*) 

AIRPLANE SPEED 0 ENGINE PRESSURE RATIO 2.04 

DIRECTIVITY ANGLE 120° ENGINE YJ75 


FIGURE 25.— F REF-FIELD SOUND PRESSURE LEVEL SPECTRA NORMALIZED TO 
1 METER RADIUS 
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NORMALIZED OVERALL SOUND PRESSURE LEVEL-dB «• » #«BAI 


REFERENCE CONDITIONS 
• DISTANCE R ■ 1 NETER 


• A r -0.0929 M 2 (1 FT 2 ) 


# p R .*.02KG/M 3 U LBM/FT 3 ) • FREE-FIELD 


•t sr -REFERENCE TENP 
(F'GURE 27 > 



DEGREES 


DATA SOURCES 
• FULL SCALE JT8D 


• HOT FLOW MODEL 
A • 43.0 CM 2 


Nj < 1.3 


440°K<T t < 810® K 


600 800 1000 1500 2000 

RELATIVE JET VELOCITY V R 
" r — — r- — t™ 


3000 (FT/SEC) 


200 


400 


600 


800 (M/SEC) 


FIGURE 26, -NORMALIZED OVERALL SOUND PRESSURE LEVEL VS. 
RELATIVE JET VELOCITY 


figure 27 -reference static temperature vs jet velocity 
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A = discharge area of the ejector 

I 

Ar = reference area = 0.0929 (1 ft 2 ) 

The 1/3 octave band SPL spectrum is calculated by adding the following corrections to the 
I OASPL, i.e., 

SPL(f,^) = 0ASPL(t) + F 2 (f/f o ) + F 3 (f/f o*^ } < 20 ' 


I 


I 


I 
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where 

F 2 (f/f 0 ) is obtained from figure 28 

F 3 (f/f Q , </>)is obtained from figure 29 

f = characteristic frequency 
o 



F^Vr,^) from figure 30. 


\/"4 

D = ejector exit diameter = VCrf) A 


When the shroud is removed, the one-dimensional flow parameters for the postmen ging noise 
region are difficult to define. This problem has been avoided through use of the one-dimensional 
flow parameters for the suppressor exhaust. This approach resulted in an empirical correction term 
being added to equation (19). The correction term is defined as a function of area ratio, AR, an 

relative velocity, i.e., 

AdB = 0.34 /AR - 1 F 5 (Vr) 

from figure 31, and the “effective” one-dimensionai fiow parameters to use in equations (19) and 
(20) above are: 


Vj = 
P = 

T S = 


velocity for the suppressor exhaust 
density for the suppressor exhaust 
static temperature for the suppressor exhaust 
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A 


discharge area of suppressor 


D = 


Vi 


r) A AR 


Premerging noise prediction. -Consider the “bare” suppressor configuration shown in figure 
23. The piemerging noire of a single element, tube or lobe. etc., ia assumed to be smubur to tha <rf a 
conventional circular jet of the same discharge area. However the mdmdual *<s for the 
multi-element suppressor interfere with each other and alter toe turbuien. s^to Cm .^dj. By 
analysis, toe affects of toe interferrring jets have been related to the number o 
elemen ts and the area ratio for the suppressor. From this analyse the space average, overall so 
pressure level for this component is defined empirically as 


OESFU - F 1 (V r , 120°) + 10 log 1Q 



+ F 6 (N) ♦ F ? (AR) 


where 

Fj(Vr, 1 2(f) is obtained from figure 26 
F 6 (N) is obtained from figure 32 
F 7 (AR) is obtained from figure 33 


V R - 


V S = 


Vj-(V s -0.2C o ) 


induced secondary velocity 

* 

flgg Cq ground static 
&|fi Cq in flight 


Cq = ambient speed of sound 


N - number of elements 

AR = area ratio: (primary + induced secondary )/primary 


( 21 ) 


. ."/* Q ; .Jf ’ ' H; . V ' ' y! iuiil'. 




DATA SOURCE: 

UAEstf PR*2^ T T *.WKimi 
ARS5.7 M,<U 



NUMBER OF NOZZLE ELEMENTS (OR TUBES) N 
FIGURE 31-EFFECT OF NUMBER OF ELEMENTS 
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and the mean one-dimensional flow parameters Vj. P, T s , A represent velocity, density, static 
temperature, and discharge area for the suppressor, respectively. 

The overall SPL varies with the directivity angle, ip , and is defined by 

GASPl (f) = OASPL + F 8 (^) (22) 

where Fg( ip) is obtained from figure 34. 

Finally, the 1/3 octave band SPL spectrum is obtained in a manner similar to that for the 
postmerging noise, except that the characteristic frequency, f Q , is typically higher and an apparent 
flttft 1 0° relative to the postmergjng noise must be added to the directivity angle. That is 

SPKf.'VO = OASPLW + F 2 ( f /f 0 ) + F 3 (f/f 0 » ^ + 10 °) ( 23 ) 

where 


F 2 (f/f 0 ) is obtained from figure 28 
F 3 (f/f 0 , >P + 10°) is obtained from figure 29 

r o ’ (^6) F 4* V R.*> F 9< M J> 

D = effective element diameter = 

F 4 (V R > ) is obtained from figure 30 

F 9 (Mj) is obtained from figure 35 

Mj = Mach number of the suppressor exhaust. 

When a hardwall shroud, L/D less than 2.5, is placed on the exhaust system, the performance 
of the configuration changes and the induced secondary Mach number increases and typically varies 
between 0.4 and 0.6 for the configuration studied in appendix B. For this short shroud, no apparent 
shielding takes place and the premerging noise level has not been observed to change significantly. 
However, this will not be the case when the shroud is lined (see sec. 5.1 .4). 
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DIRECTIVITY CORRECTIO*l~dB 




TEST DATA: 


2.54 CM DIAMETER MODEL (BSRL) 
JT8D FULL SCALE ENGINE 
J93 ENGINE 
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FIGURE 35. 


FREQUENCY SHIFT DUE TO REFRACTION AND POTENTIAL CORE CHANGE 


I 


The total noise level for the “bare” suppressor, or the ejector/suppressor conf.gurat.on .s the 
energy sum of the spectra for the premerging and postmerging noise components. The effect of 
multiple engines with this type of exhaust hardware is est.mated by use of equat.on (2 ) 

secton 5. 1.1. 6. 


5 2. 2.4 Slot Nozzle With an Augmentor Flap 


An empirical noise prediction procedure has been developed for a typical slot nozzle/ 
augmentor flap configuration shown in figure 36. The augmentor flap is unhned ** 
considered representative of that which would be used on an augmentor-wmg STOL anemft TO 
procedure is based on the data obtained from model tests (ref. 1 1 ) conducted under NASA contract 

NAS2-6344. 

Unfortunately, during the time this procedure was being developed, the recommended I nozrie 
configuration for an augmentor-wing aircraft changed (ref. 26). This new nozzle was std l a shfi, bu 
incorporated a series of screech shields in the flow. The corresponding acoustic data was not 
a vailable for analysis during the present contract period and, therefore, the application of t 
pro “l given here is limited. I. should serve as a feline for further development m no.se 

prediction procedures for an augmentor-wing aircraft. 

The current procedure is based on a static test program where the model was a 100 to 1 slot 
nozzle with an augmentor fiap positioned a, a. corresponding flap angle of — 

Table 7 gives the particular range of test conditions considered and also includes that of a 

equivalent augmentor wing configuration. 

rhe full scale equivalent acoustic data was extrapolated to the index condition (R = 1 M) and 
3 dB was subtracted to approximate free-field levels. The validity of the test data coneapon hng o 
elevation angles 0 = 70°, and 60°, was doubtful. For these data, an average between the levels 
mewiured on eac'h srte of the model was used to approximate the free-field plus 3 dB “^“ra Mo' 
overall sound pressure level data were then normalized with respect to total temp-rature ratio, 
nozzle pressure ratio, and nozzle discharge area for each value of the directivity angc 

elevation angle (0 O ). 

Tins normalization yielded a series of straight line plots with respect to the logarithms of the 

. /n T IT and area Her.ce a simple relation for the overall sound 

independent variables; P T /P S ()' T T' T TO' jnd jrtd ,Kntt * 1 „ _ ri7im T , u , 

pressure level was obtained by a least-squares-fit to the data with an R-M-S error of 1.7 dB. T 

, [/T T /P t \t / a \| (24 , 

OASPLfV', A ) ■ » 0 *’« ’°9io It-) l FTo ) K 
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TABLE 7.— RANGE OF TEST CONDITIONS 


Test Condition 

Model 

Full-Scale Eq. 
(Scale Factor*6.4)_ 

a. Nozzle Pressure Rat1o(Py/PgQ) 

1.6 to 3.0 

Same 

b. Total Temperature Ratio (Tj/Tjq) 

1.0 to 1.43 

Same 

c. Nozzle Discharge Area (A) ] 

120.6 cm 2 (18.7 in 2 ) 

0.494 M 2 (5.32 ft*) 

d. Microphone Radius (R) 

15.2 M (50 ft) 

97.5 M (320 ft) 

e. Geometric-Mean Frequency (f) 

315 Hz to 64 KHz 

50 Hz to 10 KHz 


★note* P $ q - 1.0 S. ATM. “ 2116 psfa 


T tq ^ 296 °K 


532 W R 


with the parameters: 


a 0 = function of >P and P 0 from table 8 

aj = function of 4> from table 8 

a 2 = function of 4> and 0 O from table 8 

= reference area = 0.494 (5.32 ft^) 

This relation is shown in figure 37 along with the corresponding data used to develop it. The 
norfe configuration upon which this analyses is based produced screech at noaxle pressure rafios 
exceding two. No collection was made to the OASPL to eliminate the effects of screech. However, 
ratifies i, produced in the spectra were “smoothed” in fire development of a spectrum 

shape formula. 

The 1/3 octave band, spectrum shape formula was obtained by plotting (SP1 .- OASPL) versus 
Strouhal number (ref. 7). The characteristic dimension of the slot fire, was used to calc ula e .he 
Strouha! number corresponded to the hydraulic diameter to the flow. However a temp 
stratification was observed between the hot, T T /T T0 = 1.43, and the cold, T t /Tto- ^ 
data. To further collapse the data, a modified Strouha. number was used whtch included 

temperature ratio as a factor. That is 


SPL(f.^,A>) = OASPK^A) + F(S) 


(25) 


where 

F(S) = spectrum shape curve shown in figure 38 


modified Strouhal number = f/f Q 


characteristic 


frequency = 



= hydraulic diameter “ 4 A/Perimeter 

= 2 H/(l + H/L) 


(K,L) = 


slot height and length, respectively 



















iOTES: 

/ t t V* 1 / p t V* 2 / A \ 

• $ = CORRELATION PARAMETER = 

fe) \ p so / V a r) 

WHERE ( aQ, a^ - LEAST-SQUARES-FIT CONSTANTS FROM TABLE 

• REFERENCE CONDITIONS: 

T to = 296° K (532° R) 

• FREE-FIELD, INDEX (R = 1 M) 

P so = STANDARD ATMOSPHERE (2116 psfa) 

A d = 0.494 M 2 (5.32 FT 2 ) 
l\ 
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FIGURE 37.— OVERALL SOUND PRESSURE LEVEL 
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DIMENSIONLESS FREQUENCY f/* 0 
FIGURE 38.- 1/3 OCTA VE BAND SPECTRUM SHAPE 
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Although figure 38 shows a certain amount of data scatter, the scatter is largely due t' ground 
reflection and to screech for pressure ratios greater than two. The dependence with respect to 
variations in the angles, * , and 0 o is weak. Hence, the use of a single curve is not expected to cause 
serious error in human response estimates, such as perceived noise level. 


Multiple slot/fiap conflpmtiom. -file use of multiple slot/flap configurations by STOL 
aircraft can be broken down into a single “acoustic equivalent" slot/flap configuration provided he 
spacing between the slots is small compared to the slot dimensions, i.e., less than %. e s e c 
below illustrates how this can be done. 
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“Acoustically” this is approximately equivalent to 


JM H 
T|__ 


-1 


where 


H e = MH 


L e = NL 


Thus, the use of multiple slot/flap configurations requires the change of only two input variables in 
the basic prediction defined previously, that is 


A = M N H L 


D=2MH/(I + (MH)/(NL)1 
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where H is the typical slot height, end L is the typicel slot length. 

_ _ .. „ *he recent developments mentioned in 

"ilT^rtr Nation could result through use of dilTeren, 
section 5.2.2. 1, it total tempe*ature ratio (Tt/Tjq) 

independent variables «*>» temperature ratio and eon.ec.ion Mach 

used in equations (24) and ( ) , f oredicting the noise of an aircraft in 

number. The prerent procedure has no. been .abdared f “ < **** of night could be 

flight. Taere new variables would give more insigh as obtained 

ptedicted. Further work is recommended using this approach in the analysis 

from the tests described in references 1 1 and 27. 

5.2.2.5 Externally Blown-Flap Configuration 

A mode, of an externally blown-Oap connguration is shown 

rT Howeveri tire* blown Hap has a penalty-more noise is produced 

r," a — ^ng mounting such as tha, on the cuiven. commercial airplane „ee, 

(ref. 29). 

r= z ss.’ rr: rr-'srs 

the term “impingement noise in reterence structure This results from the 

item causing an increase in noise is due to a change i the jet velocity gradients at 

presence of the flap which alteis the radiation pattern However when the 

the trailing edge of the flap can produce more none than that 

flaps are extended, the impingement noire usually “ “ the ie. noise component 

deals with the noise prediction , he “ d ditional „<** produced at the trailing edge of the 

dominates over the impingement no , refracted aft and does not contribute 

flap*, is assumed ,0 be InsigmAcan , >.<*.■£ *+ and co-annular 
to the noise observed below thc airc f . Th - exhaust as though it were a 

exhaust nozzles in sections 5.2.2. 1 and 
free jet. 

_ .. . used to de velop the nosie prediction procedure -or the blown-flap 

The acoustic f 2 8 The conditions for the static test were as listed in 

configuration were obtained .rom reterence zs. me cu 

table 9. and shown in figure 41 
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FIGURE 39. —MICROPHONE SETUP AND WING, FLAP ARRANGEMENT 


NOTES: 

0 - 5.2CH (2JD5 INCHES) 



U/D) f =7.1 , (H/0) r “13, •nr =*5 

FIGURE 40.-BASIC CONFIGURATION OF EXTERNALLY BLOWN FLAP MODEL 
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9. -TEST CONDITIONS (REFER TO FIGS. 39 AND 40) 


Item 

a. Static temperature 

b. Jet velocity 

c. Nozzle diameter 

d. Relative nozzle position 
when ® 45° 

e. Nominal Flap Angle 

f. Directivity Angle 

g. Elevation Angle 

h. Microphone Radius 

1 . Geome tr 1 c-mean- f requencl es 


Condition 


T s 

Figure 41 

V J 

D 

150 to 350 M/S 

5.2 cm (2.05 In.) 

H/D 

0.2 to 1.5 

l l /d 

-3.9 to 1.2 

L/D 

3.3 to 7.1 

*N 

0° to 45° 


10° to 180° 

fto 

0° to 90° 

R 

3.05 M (10 ft) 

f 

0.2 to 20 KHz 
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REFERENCE STATIC TEMPERATURE, T $R 



In order to extend the range of application to hot flows, a static temperature ratio term has 
been included (see sec. 5.2.2). Noise levels for flight conditions are estimated by relating the overall 
sound pressure level to the impingement velocity on the flaps and applying the Doppler correction 
for a random dipole source. Figure 42 shows the variation in the overall sound pressure level with 
the impingement velocity and flap angle. The approximate 5th to 6th power slope with respect to 
velocity implies the dominance of the dipole source for the reference configuration, i.e., L/D = 7. 1 , 
H/D= 1.5. At higher jet velocities and/or larger values of L/D and H/D, the jet noise component 
may dominate over the impingement noise. But for the reported test configuration the 
impingement noise exceeds the jet noise by 10 to 7 dB for a jet velocity range, 150 to 325 M/S 
Hence, the jet noise emitted from a STOL aircraft using the blown-flap during takeoff will 
contribute typically less than 0.5 dB to the total noise observed in the far field. 


In equation form, the overall sound pressure level at the datum condition (free-field, index: 
R * 1 M) is given by 

2. 


OASPL = [Fj(Vp, 30°) + 10 log 10 {(a * s1n 2 o^) / bj] 

1 • § / . \ 2 i 

, ) 

'sR 


where 


+ 10 Log 


10 


(Is ) V £■ ) I (1 - M cosf)- 

vt s J I'vKJ 0 5 


♦ f 3 <V'> oC n- V * W <26) 


Vp = flap impingement velocity 

* j(V r V 0 ) F 2 (L/D) + V 0 | 

F| through F 6 represent the empirical curves shown in figures 42 through 47, respectively 


V j = centerline jet velocity vector at the nozzle exit 

Yq — velocity vector of the ambient air relative to the nozzle and flap(s) 

Oxi = nominal flap angle I fig. 40) 

i 

(a, b) = empirical constants (0.0526, 1 .0526) chosen to fit the data in figure 42 


static temperature of the jet (absolute units) 
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REFERENCE CONDITIONS: 

FREE-FIELD, INDEX (R » 1 N) 

STATIC TEST, M' 0 * 0 

T $r • REF. STATIC TEMPERATURE FROM FIGURE 41. 
p R -REFERENCE DENSITY - tt.02 KG/M 3 U LBM/FT 3 ) 
Ao = FRENCE AREA » 0.0929 M 2 (1 FT 2 ) 

R 

\h » 80 ' 

P 0 -80° 

(L/D) (H/D) - 


FIGURE 42.-OVERALL SOUND PRESSURE LEVEL VERSUS IMPINGEMENT VELOCITY 




A/ *A) 


1.1 


— =6.77/ (X/D) FOR 
V i (THEORY) 


NOTES 

1) STATIC CONDITIONS 

2) FREE JET 

3) JET MACH NUMBER < 1.0 


AXIAL DISTANCE FROM JET NOZZLE, X/D 

? - CENTERLINE VELOCITY AT AXIAL DISTANCE (X) FROM NOZZLE 
j =CENTERLINE VELOCITY AT NOZZLE EXIT 
D = NOZZLE DIAMETER 


FIGURE 43.-VARIATI0N OF EXHAUST JET VELOCITY WITH 
AXIAL DISTANCE FROM NOZZLE EXIT 
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ELEVATION ANGLE, 0 O , IN DEGREES 
a) \J/ = DIRECTIVITY ANGLE RE, INLET * 0° TO 70° 

FIGURE 45 -VARIATION IN OASPL WITH ELEVATION ANGLE 







ELEVATION ANGLE,/)** IN DEGREES 
c ) \j/ = DIRECTIVITY ANGLE RE. INLET = 140° TO 180° 
FIGURE 45. -CONCL UDED 








NOMINAL FLAP ANGLE a N (DEGREES) 

NOTE 

ENGINE PLACEMENT EFFECT = Fg (X M .T) F fe (®N* 

SEE FIGURE FOR F 5 <X H ,f ) 

FIGURE 47. -F LAP ANGLE CORRECTION TO ENGINE PLACEMENT EFFECT 



t sr c 

p □ 


Z2 


A 

a r s 


reference static temperature (atraolut Jnits) trom figure 41 
density of the jet 

reference density = 16.02 KG/M 3 1 1 ibmlith 
discharge area of n© 22 le 
reference area = 0.0929 M 2 (I ft 2 ) 


M ® aircraft Mach number 
‘*0 

4 - angle between direction of aircraft motion and sound propagation path 

= directivity angle relative to engine inlet 

a s elevation angle (ftp. 4 and 39) 

K o 

X N - dimensionless engine location when ©n - 45 

= (L/D)(H/D)/8» 


D 


nozzle diameter. If the nozzle is not circular, use the hydraulic diameter, 
4A/perimeter. 


The overall SPL for N identical blown-dap configurations is estimated by adding the 
multi-engine correction, equation (2B> in seelion 5.1.1 A lo the result from equation (26) above 
The 1/3 octave band spectrum shape is shown in figure 48. Using this curve, the sound pressure level 

spectrum is defined as 


SPl(f) 


OASPL + M f /U 




where 


F 7 is an empirically derived curve (fig. 48) 

f = the characteristic Strouha! frequency in \lz 
o 

=» (Vp/D) }foMi + sin 2 o N )/(l • M 0 cos|). 
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FIGUPL 48. U3 OCTAVE -A ~ ~r r Ci P.U rf SHAPE 


The latlo, , Volo^ sin^a N ), results from Ihc change in ,he hydraulic diamc.cr of .he how 
projection on the flap(s) as the flap angle, « N . is varied. This expression assumes that the 
impingement area is an ellipse for the discharge from a circular nozzle. In reality, the impingement 
area is a hyperbolic section, and the resulting formula is considerably more complicated. However 
the scatter in the data (ref. 28), does not merit such refinement. The term, ( 1 • 0 cos{), represen 

the Doppler shift for an aircraft in flight. 

A few concluding remarks about this procedure are: 

a) The data analysis would be simplified if the reference coordinate system was relative to 
the flap(s) instead of the centerline of the nozzle, i his would permit the same data to be 

represented by feweT curves. 

b> Additional hot How tests are required to verify the methods used to extrapolate the test 
(cold flow) data (ref. 28) to that applying to real engines. 

c) For cold flow models, it is expected that the procedure presented above will predict the 
OASPL within ±3 dB provided the jet Mach number is less than one. 

d) A comparison of predicted values with a set of test data taken from reference 28 is shown 
in figures 49 and 50. 

5.2.3 Core and Turbine Noise 

The goal of reducing subsonic aircraft noise has led to the consideration of engines with lower 
„ vdodtics When compared to turbojets. The SAEje, noise procedure , ref. 7, when applied to 
these newer engines results in lower levels than that observed even when efforts are ma e . 
ehminate the fan noise component. Recent jet mode, tests (ref. 201. ,n which cure was taken to 
keep upstream noise to a minimum, have shown that a trend similar to that given s re . 
velocities less than .105 M/S (1000 fps). This observation become mme appal eij . r P 
modifications to the density correction exponent or the inclusion „l a (P T s > facto, was 
to collapse the model data (sec. 5.2.2). 

Full scale engine tests a. low power settings, however, exhibit a significantly different trend at 
let velocities below .*05 M/S than observed from the jet model tests. Apparently there are additional 
Cise sources which produce more low frequency noise than is predicted for jet 
tones due to (be last turbine stage have been identified. These sources are believed Ik gc 
upstream from the nozzle exit and can be attributed to various items: 
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PREDICTED OVERALL SPL 


t 



1) PRESSURE RATIO = 1.7, TOTAL TEMPERATURE = 283 K 

2) MEASUREMENT RADIUS = 3.05 M 

3) ELEVATION ANGLE Q 0 = 90° 

4) L / D = 7.1, H / D = 1.5 


FIGURE 49.-PREDICTED VERSUS MEASURED OVERALL 
SPL FOR BLOWN-FLAP MODEL 
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THIRD-OCTAVE BAND LEVEL IN dB RE 20 *iN/M 





Combustion 


a) 


b) Static pressure fluctuations for flow through the turbine rotors 


c) Turbulent flow impinging on the turbine rotors and stators 

d) Turbulent flow along the inner surface of the engine and nozzle walls 


e) Separated flow on the tailpipe cone and/or turbine exit struts 


Of the items mentioned, those associated with the turbine produce high frequency, broadband 
and tone noise. The low frequency broadband noise is probably due to combustion For the 
purpose of noise prediction two components will be identif.ed-(i) core noise for the low frequency 
and (2) turbine noise for the high frequency contributors. 


The present noise prediction procedures for these components are based on data provided by 
various engine manufacturers, research institutes, and NASA. Most of the data is proprietary and 
thus little substantiating data can be presented at this time to justify the procedures. In fact e 
procedures are not all that good; the tolerance is approximately +7 PNdB. A more detailed analysis 
of available engine data could result in a better core noise prediction procedure. Particular attention 
should be given to the progress made by various governmental contracts with indus ry, 

e.g., GE/FAA. 


5.2.3. 1 Core Noise Prediction 


The following prediction procedure for core noise has been developed from full scale engine 
acoustic measurements (primarily from high bypass ratio turbofans). The noise source is assumed to 
be a monopole. The strength of the source is related to the engine’s combustor and turbine inlet and 
exit parameters. Theoretically, a monopole source has a uniform omni-directional radiation pa-em 
but far-field engine data indicates that the sound is attenuated in the inlet quadrant (see fig. 51). 
The reasons for this attenuation are probably due to a combination of eftects: (l)the source is 
generated inside a duct, (2) convection, and (3) refraction. In mathematical terms, the overall sound 
pressure level for a single engine at the free-field, index condition is given by 


OASPl = 10 Log 10 


(S) 


M 0 cos |)° 


(28) 


+ F 1 (l{/) + 76 + a + K 


III 





r 


where 

m = combustor corrected mass flow 

c 

m yj T TX / T SR/( p TC/ p SR> 
m = primary mass flow 

= reference mass flow = 0.4536 KG/S (1 lbm/sec) 

T tx = combustor exit total temperature (absolute units) 

X SR = reference temperature = 288° K (5 1 8.7 °R) 

P TI - turbine inlet total pressure (absolute units) 

Pj X = turbine exit total pressure (absolute units) 

P TC = combustor total pressure (absolute units) 

Pg^ = reference pressure = 1 S. ATM. (2116 psfa) 

M = aircraft Mach number 

I = angle between direction of aircraft motion and sound propagation path 

ip = directivity angle re inlet axis 

Fj = empirical curve (fig. 51) 

a = correction for the type of burner 
= 0 for annular types 

= +9 for can types, i.e., JT8D 

K = specific engine correction (see table 10) 

The use of mulliple eneincs is lo be estimated by adding equation <2B) in section 5.I.1.6 to 
equation (28) above. 
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TABLE 10.— SPECIFIC ENGINE CORRECTION 


TYPE OF ENGINE 

jC 

CF6-6 

-1 

JT8D-9 

0 

JT8D-1 

-3.5 

JT9D-7 

+2 

RB. 211 >228 

+4 

RB. 21 1-228 

(With Revised Strut Design) 

+1 

TF34-GE-2 

-2 


Paper Engine Studies 


1 


» 


l 


l 


I 


l 


I 


I 


I 




I 


Very little data are available concerning the core noise spectrum shape. Essentially, no useful 
^ectral information can be deduced from conventional engine measurements (without a significant 
amount of work) due to the ground reflections at low frequencies. Presently, the SAE flight 
Strouhal spectrum is used, since core noise has been confused with jet noise in the past and is 
assumed broadband in nature. The 1/3 octave sound pressure level spectrum is defined as 

SPL(f) = OASPt + F 2 (f / f Q ) (29) 


where 


F 2 is the flight Strouhal spectrum shape (fig. 20b) 
f 0 = characteristic frequency in Hz 
= b/ [ V*^c ( 1 *M 0 cos5)] 

b = 1 246 Hz - (KG/sec) 0 5 = 1 850 Hz - (Ibm/sec) 0 5 

5.2.3.2 Turbine Noise Prediction 

The turbine noise prediction procedure considers two noise components: broadband and 
discrete tone. Both components have been related to the relative tip velocity of the turbine’s last 
stage, the primary mass flow, and local speed of sound at the turbine exit. The effects of 
stator/rotor spacing on the discrete tone levels is also considered. 


It has been assumed that both components have spectra shapes that normalize with respect to 
the fundamental blade passage frequency of the last stage of the turbine. The predicted spectra are 
given in terms of 1/3 octave band levels (dB re 20 m N/M 2 ) at the free-fleld. index (R= 1 M) 
condition. 


Broadband component. - The relation for the peak 1/3 octave band level at a radius of 45.7 M 
( 1 50 ft) from the source is 


L 


0 



i ®-| 0 - M cos,!) 

*"R / ° > 


-fl 


(30) 


+ F^if) ~ 10 




I 
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where 

„ Relative tip speed of last rotor of the turbine. If V TR is unknown, use 0.7 times 
the tip speed. 

= Reference velocity, 0.305 M/S (1 fps) 

= Primary mass flow 

= Reference mass flow, 0.4536 KG/S (1 Ibm/sec) 

= Speed of sound at the turbine exit. If C L is unknown, use 
C L =* a-\/lY 7 with 
a = 19.8 M/S per (°K) aS 
= 48.5 fps per (°R) 0 ' 5 


V TR 

Vr 

m 

m R 

c L 


T T7 = Turbine exit total temperature 

Cj^ - Reference speed of sound, 340.3 M/S (1116 fps) 

M = Aircraft Mach number 

o 

4 = Angle between direction of aircraft motion and sound propagation path 

if> = Directivity angle re. inlet axis 

F j = Empirical curve shown in figure 52 

sample data and predicted results am shown in figure 53. The 1/3 octave band spectrum shape 
is shown in figure 54. The sound pressure level spectrum is defined as 


SPL(f) « L 0 + F z (<7f 0 > 


(31) 


where 

f = fundamental blade passage frequency ot the last rotor stage of the turbine 
'o 

= B 0/160 (1 -M 0 cos|)l 
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DATA MAS BEEN AVERAGED AND ADJUSTED TO PEAK AT * « 110 

FIGURE 52.— TURBINE NO'SE DIRECTIVITY CORRECTIONS 



normalized SOU.® pressure level. 4,. f or turdine ROBE « » a *U 

fundamental tone 
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F 2 = Function shown in figure 54 
B = Number blades for the last rotor stage of the turbine 

0 = Shaft speed in rpm 


The use of multiple engines requires the correction , eq. (2B> in sec. 5. 1 . 1 .6, and the spectrum 
is extrapolated to a radius of one meter using 


SPL(f) 


SPL(f) 


+ 33.2 + AdB(f) | (32) 

1M 45.7 M Tabled 

Discrete tone component .- The discrete tone component of turbine noise ‘Refined in a 
manner similar to that for broadband noise. The level of the fundamental tone at 45.7 M (150 ft) 

from the source is given by 


L Si 10 L«9 10 


/ V Tr\° 6 /M / , »lVc'\(1 - M cos|) 

V~v v c l/ U/ vsJ J <33) 


+ F^(C|>) + 56 + K 


where 


C/S 

K 


stator/rotor spacing shown in figure 55 

correction for turbofans with a primary nozzle exit plane upstream from the 
secondary nozzle exit plane, i.e., the JT8D 

=» -10 dB for the JT8D 

& o dB for dual exhaust systems with co-planar exits, or turbojets 


The frequency of the fundamental tone corresponds to the blade passage frequency, f„. above. 
The higher harmonics arc assumed to fall off at a -10 dB slope as shown in figure 5<h A review °f 
available iest data indicates that the second harmonic ranges rrom 8 to 20 dB below 
fundamental. In some cases, the harmonics for the second to the last turbine stage were dominate. 
The lack of adequate information precluded further study of the phenomenon. 


In the computer program, the tones are added to the broadband spectrum (eq. (31) above) 
before the corrections for the use of multiple engines and index (R = I M) conditions are applied. 
After the corrections are made, the resulting spectrum represents the turbine noise at the free-fleld, 

index condition. 
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RELATIVE BROADBAND LEVEL- d3 











ROTOR 


STATOR (NOV) 



STATOR ROTOR SPACING EQUALS C f S 


FIGURE 55 -DEFINITION OF STATOR-ROTOR SPACING FOR TURBINES 




S.2.4 Compressor or Fan Noise 

This section deals with the prediction procedures for 1/3 octave band noise due to rotating 
compressor and fan blades. The procedures predict the noise spectra for the free-f.etd, in ex 
conditions discussed in section S.l.1.2 The methods given are applicable to turboje compressor 
single or multistage turbofans with or without inlet guide vanes, and lift fans for STOL aircraft. It i 
ZuZ that the inlets are the “fixed-inlet type," i, ; . no blow-in-doors and the blade/vane number 
ratios for the fans or compressors are optimum for minimum noise. 

The empirical noise prediction procedures discussed in this section were derived from JT9D 
and JT3D static engine tests. The data was analyzed in terms of the following noise components: 

a) Inlet fan or compressor noise emitted from the inlet duct 

1) Broadband noise 


2) Discrete tone noise 


3) Combination tone noise (buzz-saw) 
b) Fan discharge noise emitted from the fan discharge duct 


1 ) Broadband noise 


2) Discrete tone noise 

The computer program predicts each of the subcomponent (broadband, discrete tone, and 
combination tone, noise for items a) and b. above and the spectrum leve s are combinedonan 
energy basis to form a single spectrum. For an engine with more than one fan stage, each stage 
treated as an independent source and the sound energy produced by each , stage 
accordingly. No correction is made for blade row attenuation. Caution is to be applied using 
procedure for turbofans with more than two fan stages. In the case of a turbojet, the noise from 
first compressor stage is assumed to be representative of the tar-tield noise. 

5.2.4. 1 Background 

This discussion touches brietly on the various ideas and philosophies that went into the 
development of the procedures. Of the items discussed, the firs, is the definition of the 
associated with the rotating blades inside a duct with inlet and exit gun -■ vanes, c seeo 


description of some of the procedure-related corrections that are employed. Special applications 
and limitations are noted briefly. 

Broadband noise. -The term broadband noise is associated with “white” noise that is generated 
by unstable air flow past the rotors and stators of each fan stage. This noise is not strictly wh. e 
noise because at a frequency of approximately twice the fundamental blade pasmge frequency, .he 
spectral density levels have been observed to fall off at approximately 30 dB per decade. The 
corresponding 1 13 octave band spectrum shape is shown in figure 56a. The exact mechanisms which 
generate the noise are not well understood, but are thought to invo.ve local variations 

pressure field due to: 

1 ) vortex shedding from the blades and vanes, 

/ 

2) lift fluctuations resulting from approaching eddies in unstable flow, and . 

3) turbulent boundary layer(s) / 

/ 

this procedure, the broadband noise is separated into two components; one radiating 
upstream and out the inlet, and one radiating downstream and out the fan discharge nozzux Engine 
size scaling is accomplished by normalizing the inlet component with respect to the rotor diameter 
and the discharge component with respect to the exit area of the fan discharge nozzle 
normalization approximates the more classical mass Bow scaling as discussed m section 5.1 .l^The 
validity of the approximation is due to the way the turbofan engines are designed, i.e„ the inlet and 
dil^e Mach numbers do no. vary appreciably for different engines operating a. the same fan 
pressure ratio. The normalized levels are then related to the rise in pressure across each fan stage as 

shown in figures 57a and 57b. 

Discrete tone noise. - Discrete tones at integer multiples of the fundamental blade passage 
frequency are radiated from the fans and compressors of all jet engines when operated at either 
subsonic or supersonic tip speeds. A major source of the tones for fans -with in et gul e vanes 
(IGV’s) is the static pressure field developed as the blades chop through the wakes from the inlet 
JL vanes. For fans without .GVs, inlet Bow turbulence produces the same effect, bu, the sound 

produced is of a lower level. 

An additional source of nearly equal strength to the noise generated by the IGV-rotor 
interaction is the noise produced by rotor and exi.-guide-vane (EGV, interaction. A si^ificant 
parameter affecting the rotor-EGV interaction is the blade/vane number ratio. This 
Slraled in section 5.2 4.4. In both cases. IGV-rotor o, rotor-EGV interaction, .he noise is 
generated by a static pressure fleld which can be related to lift fluctuations on the blades and vanes. 
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SPECTRUM SHAPING CORRECTION 


DIMENSIONLESS FREQUENCY 
a) BROADBAND NOISE SPECTRUM SHAPE 


HARMONIC 

LEVEL 



HARMONIC NUMBER 

SKETCH NO. 1 FOR 1st FAN STAGE 
WITHOUT IGV'S 


HARMONIC 

LEVEL 


HARMONIC LEVEL = 



HARMONIC NUMBER 

SKETCH NO. 2 FOR ALL OTHER 
CASES 


b) DISCRETE TONE-RELATIVE LEVEL VS HARMONIC NUMBER 
FIGURE 56.-FAN NOISE SPECTRUM SHAPES 
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l 




FAN PRESSURE RISE.£P T /P T 
(FAN PRESSURE RATIO MINUS ONE) 


a) NOISE EMITTED FROM THE INLET DUCT 



(FAN PRESSURE RA . 10 MINUS ONE) 
b) NOISE EMITTED FROM THE FAN DISCHARGE DUCT 

FIGURE 57. —BROA DBA ND FAN NOISE (PEAK ONE-THIRD 
OCT A VE BAND SOUND PRESSURE LEVEL/ 
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In .he prediction procedure, .he dfccre.e-.one noire is separated into two noise components in 
.he same manner as .ha. done for broadband noise. The normalised levels for these components am 
shown in figures 58a and 58b. For turbofans without IGVs for the first fan stage, e g., , 

there is an Additional noise source radiating out the inlet as the relative tip Mach number ex ceeds 
!Tv Thl noise is called io.nbina.ion tones and is discussed below. As the relative tip Mach 
number exceeds one, the harmonic tones a. multiples of the fundamental blade passage 
decrease with an increase in Up speed or fan pressure ratio. This phenomenon fc though, to be the 
result of a non-linear transformation of acoustic energy from the discrete-tone noise to com ina i 
tones It is treated as follows for noise prediction purposes. At low relative tip Mach 
fundamental tone follows the solid tine shown in figure 58b. This curie is used until he Mach 
number exceeds unity After this point is reached, noise level decays for an increase in the Mac 
Zter or fan pressure ratio as shown by tile dashed tine in figure 58b. I. should be remembered, 
that the intersection point of tile solid curve and tile dashed tine vary with different engmex Th* 
point corresponds to the condition where the relative tip Mach number just exceeds unity .02 
rrrZputer progmm. The effect jus, desetibed applys only to flretfan stage without 
IGVs, otherwise the solid curve in figure 58b is used for all fan or compressor stages. 

Both the inlet and discharge fundamental tones are assumed to have their peak level at the 
blade passage frequency. The relative levels for the higher harmonics are shown m figure 56b. 

Combination tone noise (buzz-saw). -When the first rotor stage does no. have IGVs, an 
additional souree of fan noise becomes significant when the relative «P M * 

supersonic. At these hi* «P ^ * shock forms - duc , 

XZX Cwou-d be observed in the far field as a senes of tones a, multiples of the blade 
passage frequency. Experience indicates that there is a redistribution of energy. Smati diffmencc 
within the manufacturing and assembly tolerances of rotor blades appear to » « ' " 

shaoe of the shocks attached to each blade. Thus, the Mach wave system repeats itself 
revolution of the rotor, rather than with the passage of each blade. The resulting noise spec rom 
contains all harmonics of the shaft’s rotational speed. This noise has been termed combination 

noise for the subjective response it produces. 

A noise spectrum, consisting of a series of tones each with the same order of magnitude and 
separated by a fixed frequency, is referred to as a combination tone. It is a charactcns 
human auditor, system to judge the pitch of this type of noise as though it were a tone a 
separation frequency, although there may be little sound energy a. .ha.fraquencyTlvx typ 
noise ii found in all fans and compressors which operate a. supersonic tip speeds but it may 
masked by a louder tone at the blade-passage frequency or by jet noise. The presence 
mlet-guide-vanes attenuates the Mach wave system; hence, the combination tone no, sc is 
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insignificant compared to the other components as in the case of the otder by pa SS enein^ However 
if one observes a 747 aircraft during takeoff, the noise emitted from the inlets of the JT D 
is perceived to be similar to that produced by a buzz-saw. 

In the prediction procedure, the calculation of combination tone noise includes the following 
simplifying assumptions: 

1 ) Combination tone noise is emitted only from the fan inlet. 

2) It will contribute to the total fan noise only if the relative tip Mach number is greater 
than one. 

3) It will contribute to the total fan noise for fans without IGV s. 

4) Combination tone noise can be predicted in terms of three separate spec, rums based on 
peak noise levels centered at one-half, one-fourth, and one-eighth of the fundamental 
blade passage frequency of the first fan stage. 

The three peak noise levels are shown in figure 59, plotted against the relative tip Mach 
number. The spectrum shapes correspond, ng to the three peak noise levels are shown in figure 60. 


Corrections, -in the prediction procedures, various correctmns are < tm^oyed to account 

in the engine configuration from that for the reference TOD and JT9D 
change tne g. „ ro , or - s tator spacing, directivity angle, 

effects, and the use of ICVs and multiple engines. Bach 

of these corrections are discussed below. 

1 1 Rotor-Stator Spacing Correction -This correction accounts for the noise generated due to 
* Tplte TLL in front and behind of the roto. The correction is shown figure 
61 and is to be added to the peak noise levels described above in “Broadband Noise and 
••Discrete Tone Noise." Figure 62 illustrates the definition of the rotor-stator spacing 
refer to conditions at the rotor tip. 

Directivity Correction -This correction accounts for the fact that the radiation pa, tern 

for the fan or compressor noise is no, spherically symme.n. 63 nni 

procedure uses shoplifted directivity patterns for the noise subcomponents (Bg ., .63 
64). These directivity patterns show that the noise is maximum a, * - 60 fo the mle 
components and at * - 1 10° for the fan discharge components. Engine test 
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FIGURE 64 .-DISCHARGE DIRECTIVITY CHARACTERISTICS 



shown a variation in directivity angle for the maximum noise namely 50 to 70 end 90 
lo I \>° for Ike inlet and discharge components. Since commun'ty ncisc estimates are 
based on passby conditions, the erro. in EPNdB ineurred by using this simplified 
approach is small. 

3) Guide Vane Correction -The presence of inlet-guide-vanes for the first fan stage, or exit 
guide-vanes of a proceding stage, alters the observed discrete tone components for the tan 
stage being considered. The broadband component radiating out the discharge duct is ai o 
affected. The change in noise is accounted for as follows: 


a) For the inlet fan noise -subtract 6 dB from all discrete tones except the fundamental 
tone at the blade passage frequency. Combination tones (buzz-saw) are not 
considered to be significant in the far field and hence are not calculated. 

b) For the fan discharge noise-add 6dB to the fundamental tone and 3 dB to the 
broadband component. 


4) Bypass Ratio and Duct Length Comction This correction approx, mates for the change 
in the fan discharge noise that would be otaerved in the far-field if the bypass ratio and 
duct length were varied. At the present time, the physics for this effect are not 
completely understood, but are thought to be due to a change in the transmission 
coefficient at the end of the discharge duct, i.e., the duct acts as a short-wave-guide. The 
correction is defined as follows, i.e., let 


let AL 


-7.8 

6 log 10 (BPR / 10) 
0 


for BPR ^ 0.5 

for 0.5 < BPR < 

for BPR > 10 



Update AL for the change in duct length, i.e., let 

AL = (AL)C 

The constant. C equals 0, 1/6. 1/3, I respectively for short tan ducts, 3/4 length tan ducts 
long tan ducts with coplanar primary/secondary nozzle exits, and long tan ducts with a retracted 
primary nozzle, i.e., the JT8D engine. A L is added to the discrete tone levels and AL/2 is added to 

the broadband levels of the discharge tan noise. 
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5) Flight Effect Correction This it an optional correction consisting of two parts. The Erst 
part is the Doppler-shift which is well known and requires no explanation ™e 
second part includes a theoretical level correction assuming a dipole source, the result ot 
this level correction gives a slight increase in noise level in the forward quadrant and a 
slight decrease in level for ihe aft quadrant. 

An additional flight effect applies for a lift fan mounted in the wing of an aircraft. This 
effect is due to the flow distortion that results from flow separation at the inlet Up and 
the work distribution difference on the fan rotor. Since this effect is peculiar to the 
lift-fan configuration, it will be discussed further in section 5. 2.4.4. 

6) Multiple Engines Corrections-To account for multiple sources, two corrections are 
employed-one applies to the inlet fan noise components and the second applies to the 
discharge fan noise components. Equations (2A) and (2B) in section 3.1.1.6 are applied 
respectively for the inlet and discharge fan components. 


7) Index Spectra Correction-To remove atmospheric effects from the predicted spectra, the 
spectra are extrapolated inward to a radius of one meter. The far .noise pndtcban 
procedures are based on data originally measured at a radius of 45.7 M (150 ft), 
correction for spherical divergence is +33.2 dB and the atmospheric absorption correction 

is given in table 4* 

Results -Before the actual prediction procedures for compressor or fan noise are presented, it 
seems appropriate to show a comparison of the results that have been obtamcd with measured 
engine data. Figures 65 and 66 show the comparison. It is expected that the observed far-t.eld no.se 
can be predicted within ± 3 PNdB for engines similar to those in current usage. 


5 2.4.2 Inlet Fan and Compressor Noise Prediction 

This section deals with the prediction of the three noise components radiating out the inlet for 
conventional turbofan and turbojet engines. The noise from the compressor of a turbojet ^engine ts 
predicted as though it was a single-stage fan. The noise from the turbofan engme ts predtcted as the 
energy sum of that produced by each fan stage. The procedure is based on data from single an 
double stage fans. Hence, care should be exercised in using these procedures for tans with more than 
two stages because blade row attenuation is not considered. 

Hrmdbmd empouem prediction The characteristic peak 1/3 octave band sound pressure 
level of a single fan stage is defined by 
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L = F,(FPR - 1) ♦ F,(RSS) + f 3 (4>) 
0 1 & 


+ 10 Log, 


w : 


(1 - M cos£)' 


where 

^1’ ^2> ^3 represent the appropriate curves in figures 57a, 61 , and 63, respectiv y 

FPR = Fan pressure -atio, i.e., total pressure ratio across the fan 

stage being con. ered 

RSS = Rotor-Stator spacing in %, see figure 62 

* = Directivity angle re. inlet axis 

D = Fan diameter 

p, — Reference diameter = 0.305 M (1 ft) 

U R 

(1 - M 0 cos j) = Doppler-shift factor 

The result L is a free-field Level in dB re. 20 , W 2 at a radius of 45.7 M from the source at 
standard C-dihons <I5X 70% relat.ve humidity,. The sound pressure level spectrum ts 

obtained from figure 56a, i.e., 

SPL(f) • L q + F 4 (f / f Q ) (35) 


the fundamental blade passage frequency in Hz 
B 0/160 (l-M 0 cos|)l 

the number of fan blades on the stage being considered 
the wheel speed in rpm 
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Discrete tone component prediction. -The characteristic peak level in dB re 20pN/M for the 
fundamental tone of a single fan stage is given by 


L 0 “ F 1^ FPR " D + F 2^ RSS ^ + F 3 ( ^ 

♦ 10 L °9l0 [(^) 2(1 - M o COSe ' 4 ] 


where F„ to, F, represent the appmpriate ornves in figure 58b, 61, and 63, respect, vely The 
solid-line curve in figure 58b is to be used for all fan stages with the exception of the 1st stage for a 
turbofan engine without inlet-guide-vanes and operating at a fan pressure ratto greater than cnhcal. 

The calculation for this case is 

F,( F F R - 1) ■ F,(FFR 0 - l)l ol1<( , 1ne - 30-4 <-.9, 0 

where FPR ft equals the critical fan pressure ratio when the relative tip Mach number just exc 
unity The typical result of the calculation for the JT9D engine is shown as the dashed-lme in figure 
58b Tl* characteristic level, L c from equation (36), is at free-field, standard day conditions and 

45.7 M from the source. 

The next step is the accumularion of the harmonic levels to form the acoustic spectrum. The 
relative harmonic levels are shown in figure 56b. The tones am added on an energy basts *° * 
broadband spectrum. The calculation steps for accumulating the harmome levels are outtm* belouc 
The steps contain the logic for calculating only those harmonics that are necessary to tem the l/3 
octave band spectrum instead of the -brute-force" approach of calculating many 
and then determining which tones are contained in each pass imndjhe indwator GV m the lo^c 
denotes if inlet-guide-vanes are present for the first stage, t.e., IGV*0. The symbols f 0 , | 
the fundamental blade passage frequency and the cutoff frequencies for the fi. tern, 

Also, the equals sign has been generalised to denote that the results of the right replaces 

quantity to the left of the equals sign. 

Fortran instructions. — 

a o = ® * (L o ^ 

IF (IGV EQ. 0) a D - a c + 0.6 
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Nj = 1 . + f i / f 0 

DO 3 1= 1,24 

Pj = 10. **‘(0.1 * SPLj) 

N 2 = W f o 

IF ((N 2 -Nj).LT.O) GO TO 3 
DO 2 K = Nj, N 2 
IF (K.EQ. 1) GOTO! 

Pj = Pj + 10.**(a o - 0.3 * K) 


GOTO 2 

1 Pj = Pj+ 10.** (0.1 * L c ) 

2 CONTINUE 

SPLj = 10. * ALOGlO(Pj) 

3 Nj = N 2 + 1 

The results of the above calculation are twenty-four 1/3 octave band sound pressure levels for 
: inlet fan spectrum which includes the broadband and discrete-tone components. 

Combination-tone component predict (on. -Combination tone noise is calculated only for the 
st fan stage for turbofans without inlet-guide-vanes. It is assumed that only this stage ' ont " b " 

observed in the acoustic far-field. Further, flu. noise is assumed lo ears, only 

ten the relative tip Mach number is greater than one. 

The characteristic peak (1/3) octave band sound pressure levels in dB 
iquencies equal to i/2. 1/4. and 1/8 of the fundamental blade passage frequency f 0 arcg. 
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L 1 8 

a o + 

for f a f 0 /2 

4 8 

a 0 + F 5**W 

for f = f 0 /4 

4 

a o + VW 

for f = f Q /8 


where 

F 5 are the appropriate curves in figure 59 

M tr = relative tip Mach number 

a 0 " F 3 (*) 1 10 lo«li) |^-) 2 O - Mo ] 

F 3 = the directivity correction shown in figure 63 

The levels (Lj, L 2 , L 3 ) are at free-field, standard day conditions at a radius of 45.7 M from the 
source. 

The sound pressure level spectrum for this component is approximated by 

iT* inQ'Hl-K + I (38) 

SPL(f) ■ 10 Log 10 10 * * 0 j 

where G K represents the spectrum shape curves shown in figure 60 for K = i , 2, 3. 

S.2.4.3 Discharge Fan Noise Prediction 

Two components are accumulated to provide the total fan noise radiating from the fan 
discharge duct. These components are broadband and discrete tone noise. The noise produced by 
more than one fan stage is estimated by predicting the noise for each stage and the results are 
summed on an energy basis. The compressor noise contribution emitting from the engine discharge 
does not appear significant in the far-field; hence, it can be ignored. The reasons are (1) it is maske 
by more dominent sources, i.e., fan, jet, core and turbine components and (2) it is attenuated m its 
propagation through the higher compressor stages, the combustor, and turbine stages. 
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Broadband component prediction .- The characteristic peak 1/3 octave band sound pressure 
level of a single fan stage is defined as follows: 


L = F,(FPR - 1) + F 2 (RSS) + F 3 (\|» + C 

° A (39) 

+ 10 Lo9 10 C<s-> n - M o c0 4>' ,] +4 I 

R 


where 


( 


F l* f 2* f 3 ^present the appropriate curves in figures 57b, 61 , and 64, respectively. 

Fan pressure ratio, i.e., total pressure ratio across the fan 


FPR 


RSS 


stage being considered l 
Rotor-Stator spacing in %, see figure 62. 

Directivity angle re inlet axis. 

Fan discharge nozzle area. 

Reference area - 0.0929 M 2 ( 1 ft 2 ) 

Doppler-shift factor. 

3 dB for fan stages with inlet-guide-vanes. 

0 dB for the 1st fan stage without IGV’s. 

Bypass-ratio and duct length correction discussed in 
section 5.2.4. 1. 

The result, L 0 , is a free-field level in dB re 20 P N/M 2 at a radius of 45.7 M from the source at 
standard day conditons (15°C, 70% relative humidity). The sound pressure level spectrum is 
obtained in the same manner as described for the inlet broadband component. 

Discrete-tone component prediction. -The characterise peak level in dB re 20 P N/M for the 
fundamental tone of a single fan stage is given by 


A 

A R 

(1 -M 0 cos|)= 
C 

AL 
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L 0 - F,(FPR - 1) + F 2 (RSS) ♦ F 3 (*) * C 
♦ 10 Log 10 [(js-) 0 - " 0 cost)" 4 ] + Al 


(40) 


where 

F lf F 2 , F 3 represent the appropriate curves in figures 58a, 61, and 64, respectively. 

C = 6 dB for fan stages with inlet-guide vanes. 

= 0 dB for the first fan stage without IGV’s. 

The characteristic level, L„ from equation (40). is at freo-field sundard day conditions and 
45.7 M from the source. 

The nexi step is the accumulation of the harmonic levels to form the acoustic spectrum. The 
procedure for addin* the tones to the broadband component is the same as that *?***” * he 
Let discrete-tone component with the exception that we have a different charactenstic level, 0 . 

5. 2.4.4 lift Fan Noise Prediction 

The lift fan noise prediction procedure has been developed based on the fan noise prediction 

* a and a forward velocity correction 

procedure, described in previous sections S.2.4.2 and S.^,.4.3, an 

utilizing existing experimental data (ref. 30). The static noise level of an existing lift fan (ret3» 
was compared to the picdicted noise by the method outlined in the previous two sections and g 
agreement was found (fig. 67a). It is assumed that the predicted results correspond to a fan | demgn 
with an optimum Wade/vane number ratio; howerer the effects of varying the number of exit guide 
vanes or “leaning” vanes, are not considered in the prediction. 

The noise sources for lift fan propulsion systems are similar to that of conventional engines 
The fan itself constitutes the major noise source. In addition to the fan noise, none due to the jet 
and turbine are also generated by the drive system. The noise of each of these components can be 
evaluated by the methods described in sections 5.2.2 and 5.2.3. The total noise is the energy sum o 
the individual components. For lift fans that are driven by a tip turbine, tile maximum to 
Mach number of the fan will generally be less than one, mainly because of the turbine 
This limits the fan pressure ratio to below 1.3. which results m low jtf 

Therefore th* jet noise seldom contributes significantly to the perceived noise level for the lift far, 


JL 
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VELOCITY RATIO V c / V T|p 

b) LIFT FAN ROTOR EXIT DISTORTION TEST VERSUS PREDICTION 

FIGURE 67.— LIFT FAN PREDICTIONS 


a 










confutation The note doe to the ttp turbine also is of secondary importance due to the large 
number of turbine blades. i.e„ the fundamental blade passage frequency for the turbme usually bes 

above the audio frequency range. 

The effect of forward motion requires consideration of the flow distortion effect on note 
generation. This flow distortion results mainly from two sources. Flow separation rs Possible at h 
Z up due to the small radius of culture required by wing install*, ons. The second and he 
more severe distortion results from a nomuniform loading distribution existmg on the fan rotor. 
Figure 68 shows a typical takeoff flight path for a lift fan installation in two modes of operation 
Mode 1 represents vertical takeoff or the aero crossflow velocity case similar to static operation, an 
Mode 2 represents the operation with crossflow imposed on the fan. 

The complete vector diagram for the “upwind” side of the fan is shown. On this side, the 

j . a* r n Mrpcciire ratio* but on the * downwind side, the 

crossflow loads up the rotor and increases the fan pressure ratio, out on 

opposite happens. The difference in work that results can be stated as follows: 


AW~ ~ (U 2 V 2 -U, V,) 


let 


U = rotational speed - U| - U 2 

AV = change in fan exit velocity = V 2 -Vj 


Therefore 


Upwind: 


AW AV 


AW ~ - (AV + V c ) 
Sc 


Downwind: 

AW- ^(AV-V C ) 

g C 

The shore difference in work assumes that no flow entrainment takes place due to the inlet 
walls before the air enters the fan. This assumption is reasonable for shallow inlets such 


146 








I 


“fan-in-wing" propulsion system. A lift-cruise fan would have a longer inlet, and in this «, nearly 
complete entrainment of the air would take place before flow incidence on the fan face. Therefore, 
the distortion effects on noise generation could be significantly less than that cited aoove. 

The distortion for the two lift fans, for which experimental data was available in reference 31 
have been calculated. Figure 67b compares the measured and calculated distortion or the X-535 
and LF366 configurations. The calculated trend agrees well with the experiments . 
that the ratio of forward velocity to tip speed is indeed a good correlation parame 

distortion effects. 

The distortion effects on noise have been reported in reference 31 based “ f “ 

noise tests. A review of other literature and related data has indicated that the state-o -the-art 
predicting the effect of distortion is best represented by the data presented in reference 3 1 . 

The corresponding correction, to be added to the discrete fan tones, is shown in figure 63 as a 
function of the velocity ratio V C /V TIP . 

Concluding, the basic fan noise procedure described in sections 5.2.4.2 and 5.2.43 is adequate 
,0 pred r h noise for a lift fan a. static conditions. It has been shown that forward velocity ,s the 
Zl parameter producing distortion in .if. fans. Therefore, the increase in discre.e-.one noise - 
be reasonable represented as a function of V c /V T|p . Experimental dab. is currently the most 
reliable source for representing the change in noise due to distortion. 

5.2.5 Propeller, Helicopter, and Tilt Rotor Noise 

Two prediction procedures have been developed for these noise components. One is empirical 
and applies to propeller aircraft. The other has a theoretical basis and applies to hehcop.ers and d. 
rotor aircraft In the development of the latter procedure, it was found that it could also be appl 
for propeller noise prediction because the acoustic theory for propellers is esrentially the same as 

that for rotors. 

Both procedures consider two subcomponents for the observed far-f.eld noise. The 
subcomponents are < I ) discrete-, one. rotational noise and (« broadband vortex > 
procedure the vortex noise is predicted by empirical equations because the more refined integration 
and boundary-value problem approaches are computationally expensive and '^ ^Ter/Tr 
information than is readily available. These refined approaches are imporlant to proMc / 
design; but the increase in accuracy for uWurc levels is no, .ha, impressive. The two procedures 
described here differ in only one respect. The rotational noise is predicted empirically 

theoretically in the other. 
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An error analysis has been made comparing ,he ,wo procedures with daia^.e^in recces 

s — :rr: 

for the tea, conditions with the firs, destructive interference occumng between 

The data analysis is too lengthy to present here.hu, i, is worth noting the results that were 
. KiarirsH H R 212-14 propeller. Based on a sample of 90 spect.a, the 907 
obtained tor a f °J; 7 . 0 y B and (+5.4\ pNdfi for the rotor and propeller prediction procedures, 

study are required. Some sample predictions are shown in figures 70 and . 

5 2 5 1 Propeller Noise Prediction 

Sr—SSS-SS r#r 

meln't references 33 through 35t bu, i, is rarely present for conventional propelier atrcraf, 
operating at subsonic tip speeds. 

The procedure defined for vortex noise is based on a combination of the empirical approaches 
P n tiv rti.oh 3R Reference 32 gives a simplified, empirical procedure for 

Vortex noise. The equation for tne overall sound pressure level, L Q , in dB re 20 pN/M at a 
distance of 1 52.4 M is given by reference 35 as 


20 log 10 [^ V TE T 1 ^ V R V 


10 log 10 


* A R #0.1 ", 

K W 


cos 


0.1 + cos 70' 


It) 


(1 


M cosl 

0 
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FREQUENCY IN HERTZ 


CONDITIONS: 

• FREE-FIELO, INDEX ( R = 1 M > 

• TIP MACH NO. = 0.4, SHAFT SPEED - l7 j° “ $pAN 

• FOUR 5 FT. DIA BLADES. BLADE ANGLE - lb AT 0 _” SP 

• THRUST = 324 LBF, TORQUE = 147 FT-LBF, SHAFT POWER 


48.1 HP 


FIGURE 70.— MEASURED DATA IH.S.212- 14 PROPELLER! VERSUS PREDICT, ONS 
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pww MEASURED 1 DATA SCATTER FOR GROUND MEASUREMENTS MADE AT 
U “* J | VARIOUS ANGULAR POSITIONS ON THE 61 METER POLAR ARC 

pm PREDICTED J 



TEST CONDITIONS 
ISA CONDITION 
HOVER 

AIRCRAFT WEIGHT • 20,000 N <4,500 LBF) 
BOELKOW-105 HELICOPTER 


FIGURE 71. -sample helicopter prediction 


V 


* 


UL 


where 


V TE 


= effective helical tip speed 


Vr 

My 

M„ 


C 0 -^M T 2 + (M 0 cosa) 2 

reference velocity = 0.305 M/S (1 fps) 


= tip Mach number 


a = 


aircraft Macn number 

angle between propeller axis and direction of aircraft motion 


r = ambient speed of sound 


T = tluust developed by propeller 


reference thrust = 4.45 N (1 lbf) 
total blade area on one side of the propeller 
reference area = 0.0929 M^ ( 1 ft^) 

* = directivity angle re inlet axis of propeller 

5 = angle between flight path and sound propagation path 


Tr = 
A b * 

a r = 


However, both references 33 and 35 point out that this fonnula pertarns to prope lens wrih .5 
or 6 blades and a correction of *5 dB should be added for conventional propellers of 2 to 4 bl ades. 
■Mso corrects ms for atmospheric absorption and (pound rellecUon appear not to have been me 
References 33, 35, and 38. I. is desired to have the overall a. free-field index condritons (R- 
1 M). Hence, the constant -43 dB in the above equation is adjusted as follows. 


-43.0 

+42.7 

+3.0 


from original equation 

spherical divergence between R| = 152.4 M and R 2 - 1 M 
estimate of atmospheric absorption at 1 52.4 M 
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Y 


v 


+5.0 


number of blade correction 


-3.0 


correction to free-field condition 


+5.7 Total 

The resulting equation for the overall SPLofa single propeller at free-field index, conditions t> 


L 0 - 2C 1.5g 10 L* TE T / <V R T„)] 


* ’0 ( 0 ~! * 0 - M o **>'''* 5 - 7 <4U 

10 m B 0.1 + cos 70° 


The characteristic Strouhal frequency f 0 is given by reference 37 as 
f = 0.28 V/l(t cos 6 + 4 sin J) (1 - M 0 cos$)l 

where (V t « , d ) are the effective helical 'velocity, blade thickness, chord length, and angle of 
attack respectively at 0.7 span. The basic spectrum shape, SS, which includes frequency modulatron 
effects of blade rotation is given by the formula: 


SS * 10 tog 1Q 


h 


[< 2 «(*) <1* 

k 2 


] 


r [1 * (a S 2 ) 2 ] [1 ♦ <b 
|[1 ~ (b S./] [1 * (» s,)^! 


(42) 


where 

S] “ fi/f 0 

§2 ~ ^2^o 

f , , f 2 = lower and upper cutoff frequencies for the pass band being considered 


g(x) = 


the normalized, power spectral density function 


f_aL_. _b L_“ 

1(1 + a 2 x 2 ) (l+b 2 x 2 ) 


1 +M 


TE 


b = 1 * M T E 

c = log e (a/b) 

j^Te = effective helical tip Mach number. 

The spectrum shape formula has a singularity at M TE equal to one and thus the pr ° c ® d “ r ®_J^ 
. J2E «P speed soes supersonic, vortex noise ^icioo »' " ' ^ 

bands to the overall sound pressure level, L Q . 


^Pl (f\ - l. 


+ ss 




Rotational noise. -The characteristic level in dB re 
noise is given by reference 32 as 


20 p N/M 2 at 152.4 M for the rotational 


L 


o 


«r ]•■* 

2.2 B + F 1 Cl|0 + 100. 


where 

W = shaft power 

Wo = reference power = 745.7 KW (1000 Hp) 


O - propeller diameter 


= reference diameter - 0.305 M (1 ft) 

B = number of blades 

F| = directivity correction (fig. 72) 

The expression shove is a — to the 

that this level is representative of measure groun reflection would be between 2 and 

frequency for destructive interference ue o . j spec trum. This implies that the 

5 KHz-well above the first ten harmonics in the rotational ^^ s 33 and 36 report that a 
... ii ie « HR above free-field conditions. References 55 ana oo rc F 
charactenstic level L 0 is 6 d , tiftna i motion. Thus to meet the 

— ^ characteristic level a, «he free-f.eld. index 

condition. r , , lV 1.55 /* \ -2.265 


10 Log 


38 M te - 2.2 B + ♦ 137.7 


0 


- M q cos£f 4 j 


(44) 


/i m is the Doppler factor. This last expression assumes that the air 

" "e^le due to J fundamental frequency for the propeller tones being 

typically less than 250 Hz. 

The levels of the hamtonic tones are detained thrornh use office 73 for the function F 2 
used in the equation below. 

.. l 0 ♦ K). K - 1.2.3. ••• (45) 

The fundamental frequency for the first harmonic is defined as 


f . B 9 / [60 (1 - M C0S&] 
o 


(46) 


Where i Is the shaft speed in rpm. The discrete tones am added* the broadband spectrum in a 
manner similar to that described for the fan noise procedures (sec. S. ..4). 

5.2.S.2 Rotor Noise Prediction 

A — d medic, ion 

procedure considers two noise components that are general 
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I 


HARMONIC NUMBER 


are ( I > discrete-, one, rotational noise and (2) broadband, vortex noise. Rota tonal no.se ts 
developed from the harmonie loads that exist on the blades and that for vortex no, se ts due to the 
shedding of vortices from the trailing edge of the blades. Other noise sources are ment.oned 
references, but they are usually neglected. The reasons are discussed at the end of thts sectton. 

Vor,e X noise. -The previous section 5.2.5.1 provides a simplified, empirical procedure for 
estimating this component. The method used here is identical to that descnbed tn that sectton. An 
ereor analysis was made on some of the data given in references 37 and 38. The results tndtcate a 
the empirical method provides reasonable predictions for communtty notse eshmates and it 
“cheap" in comparison to refined integration and boundary-value problem approaches. Also, they 
require more information than is readily available, e.g., the complete blade geometry, etc. 

Rotational noise. -The procedure for rotational noise is based on a theoretical math model 
(ref 33) that is stmplified by a "loading-law" concept (refs. 35, 37, and 38). For the purpose o 
noise calculations, the harmonic loads are considered random in phase and apphed at a stngle prnnt 
on a blade. The position of tins equivalent point load is based on a centroid calculahon f 
mean-square pressure distribution on the blades. Experimental data (refs. 35, 37, an ) 
shown that the load harmonics can be estimated from a simple formula. This formula contains on y 
three empirical parameters which can be determined from acoustic data. 

Reference 35 gives data which show that the harmonic loads are different for a rotor wtth 
translational motion Hence, the loading parameters will have to i* obtaine torn wtn tunne 
tests in order to predict the rotational noise for Right conditions. From what little data t 
available, the discrete-tones for flight could be significantly less than that from a stahc atrcraft 

(ref. 37). 

One may ask. "Why don’t you assume some form of a time-varying pressure distribution on a 
blade; compute the Fourier series; and use the coefficients thus obtained for the loading 
harmonics?” The answer is-it was tried and wasn’t successful! (ref. 33). Furthermore,! ts t ta : er 
computationally expensive and each rotor design requires a different pressure distnbuttom 
On-the-other-hand, the loading-law concept simplifies a rather formidable ma,heraa *“l^ ^ are 
and provides realistic noise estimates. At the same time, it lumps together many rife* that are 
difficult to predict into oniy three empirical parameters. Some of the effects me blade fiapp.^ 
flight, vortex interaction, and changes in blade design. Of the various means avatbble to measure the 
harmonic loads, the acoustic method (refs. 37 and 38) seems to be the best. This ts because othc 
techniques (ref. 33) are presently unable to determine the harmonic loads to the htgh-order requtred 

for noise prediction. 
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Rotational noise for helicopters and tUt rotors is composed of discrete-tones that occur at 
harmonic multiples of the fundamental blade passage frequency, f 0 . 


f Q = B 0/(60 SF) in Hz 
B = Number of rotor blades 
0 = Rotational speed in rpm 

SF = Doppler-shift factor (1 - M c cos£) 


M 0 = Aircraft Mach number 

| = Angle between flight path and line to the observer at “retarded time”-the time the 

sound is generated, not the time when the sound is heard. 


The harmonic levels, dB re 20 PN/M 2 for the rotational noise are given by 

2 
C, 


L n = 10 Log 


10 


N 

^Tp, 


= 124.57 + 10 Log 


10 


with 


4> 0 = 20 PN/M 2 = 4. 1 77 x 10' 7 psf 

P R = 47.88 N/M 2 = 1 psf 

Reference 33 provides a theoretical approach for estimating the far-field, acoustic pressure 
produced by harmonic loads on a rotor. The result, equation (34) in reference 33, gives an 
expression for the discrete tone phasor, C N , above in terms of these loads. 


C N 3 C TN + C DN + C RN 


(47) 


I 




I 


I 


i 


* 


I 


4 


4 
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'TN 


A .-(n-A) / infix \ . , /jiMy 

•£.' Weir) TX \ r i 


• Thrust component 


^ . -(n-A) / i(n- A) \ / nMi 

c on = \2 wrjR / A ° x r, “ x V r x 


nM y ^ 


a Drag or torque component 


C,. . £ i 

■ A .-» 


s. .-(n-X)/ -nv \ . /^i\ 

\2*c,r|j R * "" l \r, / 


Radial component 


where 

n = NB 

rj = rSF 

r * Distance from observer to rotor hub at retarded time 

M = Rotational Mach number of a point on the rotor 

= 0R/c o 

q = Rotational speed (rad/sec) 

R = Radius of a point on the rotor 
c = Local speed of sound 


* = COS Ip 

r 


L = snip 


ip = 


Directivity angle between the rotor inlet axis and a line to the observer at 
retarded time 


i = Imaginary number, V-T 

(A Tx . A dx , An) the complex Fourier coefficients for the harmonic loads appearing 

on a rotor blade, i.e., 

. . JL f f (T) exp (-IXQT) dT 

X 2TT •£ tf /Q 


f( T ) = ^ Ax exp (iXQT) 

However, reference 35 point out that an “effective” helical Mach number, Mg, is to be used 
in lieu of M above to account for motion relative to a stationary observer. 

Mg 2 — M 2 + (M 0 cos a ) 2 

in this expression, a is the angle between the rotor axis and the direction of the aireraffs 


'Th 


'DM 


C R'I 


t Vxl 

2URrSF & 

( nM E sin \K 

^ SF / 

-InMg A , 

S OX 1 Mp \ n 

2l?RrSF £ L E 


nMpSln ^ 1 A J* 

— ; L 1 RX n=A 

2 lTRr$F z 

/nM E sin^\ 

\ SF / 
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Changing the order of summation from (-«*•) to (1,-) and noting that A„ x = A x *. complex 
conjugate, yields 


C N = C TN + C DN + C RN 


(48) 


with 


1M e cos<|» 
” 2itRrSF 2 


[■ 


a T0 


09 -(n-A) 

♦£ ■ | A TX J n-A * ( 

A«=l 


- 1,A A n J n*x] } 


-InM* 


'DN 2wRrSF 2 






r 2 
l RN 


A 

nM E slni|» 


2^Rr SF 2 


f ! - n a / 
— fl A Ro n 

SF 2 l 


00 -(n-A) r . » . » 11 

♦ £ i I A RX I n-A* ( " 1) A RX 3 n^x || 

An L J 


where the argument in the Bessel functions above is (n M E sin*/SF) and the following identities are 
to be used in evaluating the negative order Bessel functions and their derivatives. 

^n-A ” [^n-A- l^n-A* J 

J n4A =0oS [^-t-A-l^n+A+l] 


It 


Loading laws .- Up to this point, it has been assumed that the harmonic loads (A JX , A DX , 
Arx) are known; in reality they are not. References 33, 35, and 36 argue that the harmonic loads 
are to be considered random in phase with respect to harmonic order, X, and position, R, across the 
rotor. For the purpose of noise calculations, an equivalent load is assumed to be applied somewhere 
between 0.5 and 0.8 span. The position is based on a centroid calculation for the mean-square- 

pressure distribution on a blade. 

This point load concept permits a simplification of equation (48), i.e., the values of M E and R 
at the centroid can be used, to relate noise to the steady forces on a rotor. The harmonic loads are 
in turn related to the steady force by the following approximate loading law. 

2 |A X | / A q = g(\) (49) 


where LI-(m+0.5) 

g(x) - c |*| 


and m and c are determined by physical argument and/or experiment. 

These loading laws imply that an effective phase for the Fourier coefficients, A x , is 45° for 
X>0. Also, the summation in equation (48) with respect to X is to be done on an R-M-S basis after 
the algebraic sum of the load components. Free-Held, acoustic measurements on the inlet axis of the 
rotor, theoretically, provide an estimate for | A TX | with <r > 10 max[c 0 /(Nf 0 ), Rjl ard X = N B). 


That is: 


M 2 


it R r SF 2 p Q 


XMi 



Similarly, measurements in the plane of the rotor can provide an estimate of |A dx | 2 if the radial 
forces are small in comparison to the drag forces due to torque-i.e., 


A RX | <<| A DX 


See references 37 and 38 for furthe 'tails. 

References 33 and 37 say that c - 1 and ~ 2 for hovering helicopters. Thus 


, , 3 (50A) 

9(A) - M 
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On the other hand, reference 35 gives a physical argument and data (fig. 74) which show that 
equation* (49) applies with c = 1 and the exponent term, m, determined by a « winch 
includes effects of rotor orientation and aircraft speed. Theref- for helicopters and fit rotors. 

, -(m+0.5) 

«w-M (50B) 


where 


~ 0.0485 


ri6ir(V 0 co s >or+v)1 , , 

l aS SI Rj J “ 


= 1.3 + 0.4876 IVjS^CSaVj)] 


V 0 = 


Aircraft velocity 


Induced velocity of the air 


~ 0.5 |*(V 0 cos ) + l V Q cosoc) 2 +<i2 g c T)/(p 0 % 2 )1 } 


V, = 


Normal inflow velocity 


/ 2 T so T *"» \ 

« 0.5|(V o cos )+[(V 0 cos ) z + a — - 2 J } 


so u T 


Vj = Tip speed 

T = Total gross thrust for rotor 


Density of the air 


Pso/< R cT S o> 


Static pressure of ambient air (absolute units) 
Static temperature of ambient air (absolute units) 


Tip radius 


Tip diameter 


it 






I 


i 


tj = Blade area/disc area 

S = Lift-curve slope of one rotor blade 

S = Reference lift-curve slope - 3 for data in figure 74 

o 

a = 730.97 (M/S) 2 per (*K) 


I 


f 


I 


I 


= 7.2141x1 (T 3 (M/S) 2 per(°K) 


N/M^ 


= 4371.03 (fps) 2 per (°R) 


Be = 


i.oS!l“ 32n40S ^T 

N-sec^ lbf-sec^ 


R, 


M7 05 j^= 1833* lOrJ^r- 


- 53 35045 kpr 


For low-speed, Mj 


< 0.3, propeller noise estimates, reference 37 gives the empirical loading 


law 


-1.43 


g(X) = 0.86 |x| 


(500 


and reference 36 says; for the same tip speed range 

g (X) * 0.04 [X 2 (1 ♦l x /36| 3 )] 0,5 

There seems to he a typographical error here, because this fommia doesn’t match the empirical data 
in reference 37. The correct formula should be 

, -I -0.5 

(SOD) 

.■ i th ■*' i i .iir j w ■ 

9(M 


1.22 | 1+(X/36 * j] 

44.0 |xr 2 V[l4(36/X^] 


0.5 

0.5 


I 
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Although the two equations, (50C) and (50D), appear quite different, they match the same 
data (fig. 75a) and have one thing in common when compared to the equations for helicopters- the 
propeller harmonic loads are typically higher. 

More recent data (ref. 38) at the propeller conditions, 0.5 < M T < 0.7, yield yet another 
loading law 

g(X) =s 24.4 N' 2 ' 5 / Cl ♦ (30/X) 2 ] 0 ’ 5 (50E) 


This formula is derived from the data given in reference 38. The spread in the data between that for 
the thrust and torque harmonics is probably due to-the equivalent point loads for thrust and drag 
act at different centroids. Thus, the use of a single point load for calculating the harmonic tones 
results in a tolerance of about ±8.0 dB. A plot of equation (50E) is shown in figure 75b. In view of 
the data and formulae presented, one would expect that the general form of the loading law is 

n -0.5 

(51) 


with the parameters (c, r i,X c ) determined by experiment. In actual practice, the form of equation 
(49) should suffice for determining the sound harmonics of orders; 1 < N < 30/B. For more 
accurate estimation of sound harmonics outside this range, equation (51) should be used. 

Simplification. -In order to put equation (48) in more useable form, the following notation 
is used. 


a to 

= T/B 


a do 

= h D A to 


a ro 

55 h R A TO 


| a tM 

= 0.5 A to g(A) 

for |*| > 0 

| a dx| 

= h d| a tx| 

for |X| > 0 

| a rx| 

f< 

< 

II 

for|x| >0 
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i-propeller harmonic loads 






«. . — «» - » - - - «“ “2T,r“ ST” *"* 

(50E), and (51 ). Assuming random phase, as was done in references 33 ar. , 

A k = !ii! ( 1 ♦ 1). X> 0 
y/Z 

for each load component. Insertion into the random phase form of equation 48 yields 

' -• ** ^ \2 r 


N* =(S)l 2( “ 


2 +8 2 ) 
no **no ' 


♦ T^V) (a nX 
X*1 


'nX 


] 


where 


1) X Jn+X 


“nx = [ h R $ln * J »-x] + [ c45 *‘ l ’0 («e) ( n 

*.X* [ cos* -« D Jn-x- [»« s,n * <-*> X 
(n M e sin*/SF) = argument in the Bessel functions 


n 


NB 


M e 


Helical Mach number at radius R 
Radial centroid for equivalent point load. 


(52) 


Equation (52) abo« provides an estimate of the MU discrete tones a. a distance, r. from 
the rotor huh. This estimate also corresponds to free-fleld conditions. 

!dd to the noise. References 33, 37, and 38 show that the effective rante for X B 






with 


q <= |m e sin*/ SF | 
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and the terms containing J n+ x»“<l 4+ X can ** neelected ' 


„ also can be shown that the radiai load components can be neglected in equation 52 because 
their amplitude is much less than that due to thrust or torque. This results m a less complex 

equation 

_ . I f a v V- -> 

Q) 




2^9 2 <x> J n ?x< n( i> | c °s*-ho(^)( iL jr)j J 


where 


*1 - 

K 2 = 

q = 


Max(n (1 - q) -0.5, 1 ] (integer result) 


n (1 + q) + 0.5 


(integer result) 


Mg simt/SF 


n 


= NB 


Mg = 
R = 
g(A) = 
T = 


Helical Mach number at radius R 


Radial centroid for equivalent point load 

Represents the loading-law function, equations (49) through (51) 


Total thrust for rotor 


h D = Drag/thrust ratio 


$ = 

SF = 


(Q/R)/T where Q is the total torque on the rotor 
Distance front rotor hub at retarded time 
Directivity angle re. rotor inlet axis at retarded time 


Doppler-shift factor, 1 - M Q cos | 




(53) 
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The use of equation (53) to calculate all the sound harmonics required in the audio fluency 
.JZZ be Z* expensive-even on h*h«peed con.pu.ers. Reference 33 has shown toa. .he 

«u!tine values of L w vary in a smooth fashion when plotted against log(N). Suppose that a data 
“ 2Z to can be developed through use of equation (53, for a selected set of 

values for N. Then by means of interpolatton/ex.rapola.ion wilh respect <o log(N), M the necesmty 
I" harmonic I JL I* can be obuined a. a tendons saving in computer storage and t,me. 

Additional limitations must be employed to order to .educe computer storage and tone. The 
practical limitations for rotors and propellers are 

2<B <6 
Mg 

0< SF <1 

As was mentioned previously, the values of N will be limited to 

1 <N <21 

Also, symmetry implies that the values for sin* are contained in 

0 < sin * < 1 

Thus the maximum range for X in equation (53) is given by 

1< X <2n 

where n = N B and the orden. of the Bessel functions, J K , which could have to be calculated and 
stored are 

0 <K <n 


Inserting the possible values give 

0<K <42 for B = 2 
0 < K <63 for B = 3 

0<K < 126 for B = 6 
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RELATIVE DISCRETE-TONE LEVEL (l^ - L t ) 


l 


i) ' 




or in general, 


0<K< NB 


with N being the harmonic number for the discrete-tone occurring at the frequency, N f 0 . 

Other harmonic noise sources-There are three other noise sources resent for gating 
and/or rotors. These sources are (I) "thickness-noise” (dipole character), (2) Reynolds 
Stress noise (quadruple character), and (3) blade-slap. In practice, these no.se sources are usda y 
neglected. The reasons are discussed below. 

In addition to the thrust and drag loading noise already discussed, another dipole source is 
present (refs. 36 through 38). This source has been given the term “thickness noise dm to _.ts 
Length being proportional to the blade volume and the local blade acceleration (re ). 
rotora operating at constant speed ( 6 and M 0 constant), the only acceleration present^* m he 
radial direction. (Note that the blades am considered rigid.) Thus thra srmrce 
radial loads that appear on the blades. In order to evaluate this source, knowtedge of the comphde 
Made geometry is required to perform die necessary integration over the 
Standard reports (refs. 37 and 38) have included this source in ‘^ “'^''“^J^gh.ly 
this source can dominate over the loading noise due to thrust and drag if the 
loaded or if the blades are rather thick. Reference 36 argues that #us souree wdl only ^be dgm • 

relative to the thrust and drag ternis, for very thick blades at low loading conditions^ Although th - 
three reports agree in concept; the emphasis of reference 36 differs. This latter emphasu; leads to the 
assumption: “The thickness-noise contribution to the discrete-tone levels produced by propellers/ 
rotors can be neglected, if the blades are not thick and/or lightly loaded. 

The relative magnitude of the quadrupole sources when compared to that for dipole sources is 
given by reference 36 as 

(c„ / C N I) - »"($)«»« 

Quadru- Dipole 

pole 


This relation shows that the dipole (force) noise components dominate in most practical 
applications. However, as M cos 4 approaches unity and n gets greater than 100, the quadru pole 
sources become dominant. “In assessing these results, it is of course nnportant to re^mto ^ 
these results apply only to the specific case examined. . . .Despite these irniitatrens .. may be 
concluded that for quiet propellers operating at tip speeds of less than about ” ' ' ’ 
quadrupole noise should not be expected to make significant contributions to the harmonic 


for n < 200, even on the propeller axis” (ref. 36). For further details see pages I, 9, 68 ot 
reference 36. 

Under the various helicopter operations (for instance, during low-power descent), the rotor 
produces a loud impulsive noise. The energy for this noise consists of harmonic tones that occur at 
multiples of the fundamental blade-passage frequency, but the distribution does not fall-off rapidly 
with increasing harmonic number-hence, the impulsive character. This noise has been given the 
term “blade-slap.” Whenever this “impulsive” noise occurs, it is particularly severe, but it also has a 
highly directional radiation pattern. Thus blade-slap is m i always heard, even though it may exist 

(ref. 33). 

This phenomena occurs at precisely those conditions when a vortex wake can be expected to 
pass very close to the rotor. Blade-slap can also occur when the rotors are operated at high speed. 
Then, it is associated with transonic flow over the rotor blades. Thus, there are two possible sources 
of blade-slap-vortex interaction and/or transonic flow. These phenomena can be predicted. 
“However, it seem inappropriate to consider blade-slap as a separate phenomenon. The helicopter 
rotor is always undergoing some form of vortex interaction, and blade-slap is simply a severe form. 
Perhaps it is more realistic to suppose that, at least from the acoustic point of view, the helicopter is 
always flying under some degree of blade-slap” (ref. 33). 

The loading-law, equation (Sot), includes some of these vortex interaction effects and its use 
is recommended. Hopefully, this equation approximates the low-degree, blade-slap mentioned 
above. No additional effort was made, however, to try and predict the special severe case. 


S.3 NOISE CONTOUR ESTIMATION 

A noise contour is the locus of points on the ground in which the noise is at a constant 
acoustic level. The calculation of a noise contour requires the establishment of the relationship 
between the aircraft’s noise performance and the aero/propulsion parameters during takeoff and 
landing. The following optimized method is presented which will fit within the computer time and 
storage constraints of the Ames flight simulirtor. 

The relationships mentioned above are established when data points are given for noise level 
(NL), engine performance (EPP), range at closest point of approach (R), and elevation angle (a), for 
an aircraft during level flight (see fig. 1 5 and table 11). This data can then be formed into tabular 
functions: NL versus (EPP, log R.or ), or log R versus (a. EPP) for each noise contour. When the 
airplane coordinates and EPP are given, interpolation using these functions at the geometry shown 
in figure 77 provides two points (one for each side of the flight track) on the ground for a specific 









FIGURE 77.— NOISE CONTOUR GEOMETRICAL RELATIONSHIPS 


noise level. .1 a series of these points a. catenated during an aircraft's takeoff or landing, a noise 
contour is determined. The area enclosed by this contour can be calo.la.ed. 

Although there exist more refined methods fo, calculating noise contours, they req uire rater 
tengthy calculahons. and Z has the advantage of 

minimizing computations and reducing storage requirei^irts^Des^e^^fact ^by^fned^ction^the 
approximations, when sufficient data points are provided by measuremen 
procedure provides reasonably accurate noise contours. 

5.3.1 Acoustic Data 

me acoustic data required for noise contour estimation consists of a dtotivhy angleta ** 

noise radiadon and a tabulated function of three variables-no.se M < 

fie 15) It is of particular importance, when constructing this function, that file data a tor e 

source estimation procedure. 

i- SJ.2 Aero/Propulsion Data 

The aero-propulsion data required for noise contour estimation consists of a series of points 

a, ,he " ^ 

orientation angle (6 E > tor me ruere 4 . , „ ntMCt the key engine performance 

.. m m , FP m Durine the Phase A portion of the contract, tnc Key uiy v 

performance (LP ). D g relationship to other jc. engine cycle 

parameter was engine pressure ratio (EPR) d 

* , : „ thort. is one-to-one correspondence between EPR and all otner engine cy 

parameters and this »» this 

no.se pioduced by jet engines is a.re y However ict engines are not the only 

”"e d' n thl Phase B effort. The choice of wha, the engine performance parameter 

represents is left up to the usei . 


5.3.3 Noise Contour Calculation 


The noise contour calculation can 
formation of the acoustic data functions: 


be broken down into four basic steps. The first is the 


and 


NL= fj(EPP, log R,a) 
Log R - $ 2 ^*' 


for each noise contour from the given data points (NL k , EPP k , R k , « k >. This step ts done on y once 
Thereafter, the caicuia.ion requires interpolation for iog R a, a desired contour no.se ievei (CNL) 
"letry shown in figure 77. The next steps are the geometric soiutions for fc — p« 

(U V) in a moving reference frame and finally, the transformation of contour coordinates (U, V) t 
the fixed (X, Y, Z) coordinate system (see fig. 77). The details are outlined below. 

a) Formation of the acoustic data functions, fj and 

Dam points (NL k , EPP k , R k ,« k > are assumed to he given " - of the noise 
source estimation computer program or from measurements (table 1 1 ). Sort the given 
data with .aspect to increasing values of EPP k , R k , o k as these variables will be treated 
independent for the function f,. Next, determine the distinct values for the given arrays 
(EPP k , R k ,o k ) and use the results for the independent variable data arrays speci ying 

where noise levels are defined. 

NOTE: The three-dimensional function NL = f j(EPP, log R,«) is now formed. 

Specifying the desired contour noise levels (CNLj) permits the transformation of the 
function f, to a function log R = f 2 («. EPP) for each CNLj. This is done by 
one-dimensional interpolation on fj at NL - CNLj for j = 1, 2, etc. 

NOTE: The transformation assumes that the function f , is monotonic: decreasing with 

respect to increasing values of log R 

b) Calculation of log R for a specific contour CNLj: 

At each point along the aircraft’s flight path the folio ving data is required. 

Zj aircraft height above the ground 

6 . climb angle, i computed as the are tangent of the climb gradient 


6 Ei 


Orientation angle of the dominant noise soutce reference axis-usually the 
angle between the gross thrust vector and the horizon 


EPPj engine performance parameter 


Calculation: 


let Z E = aircraft-to-ground distance perpendicular to the flight path 


Zj COSdj 


Iterate the calculations below until |a-a 0 |<(c o|. In this iteration, a 0 is set initially to 
r:,le 45* for each no.se contour and'is Updated at each aircraft pos.hon along the 
flight path. The value, e, is a tolerance constant for the iteration; a reasonable value is 

1.2 x 10' 3 . 

log R = interpolation on f 2 at (o 0 , EPPj) 

Of - arcsin (Z E /R) 

Test a for convergence with and update 0f o if another iteration is required. 

In the computer program, a test for contour closure is made just before the calcU “°^ 
a above. Closure occurs when R < Z E . The action taken is to set a - 0.5 + 90 ). Tta 

is done to avoid premature closure estimates that have occurred during jmij b 
operations. If closure has indeed occurred, the program contains a trap 
corresponding error code. 

c) Calculation of contour points (U, V): 

U 2 = r 2 -z e 2 


V = Solution of 


I 


P cos $ = v cos d Ei -Zj sind E j 


' cos d Ei -Zi sind E i \ 

U 2 + V 2 + Zj 2 J 


where 


R = 
* = 


range at CPA calculated for the noise contour in step (b). 

directivity angle for peak passby noise propagation relative to the 

dominant noise source reference axis. 
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NOTE: Singularities can exist when fig, = -90 , and/or the directivity cone does not 

intersect the ground at sideline distances ±U. Helicopters, tilt rotor aircraft, etc. 
require special consideration. The singularities can be avoided by letting 

dfij = and ♦ = 90 . 

d) Coordinate transformations: 
let dx = x r XM 

dy = Yi-yi-l 

dsj 2 = dx 2 + dy 2 

sin 0 = -dx/dsj 

cos 9 = dy/dsj 

x = U cos 0 - V sin 0 + Xj 

y = V cos 0 + U sin 0 + yj 



for each contour CNLj where (from section 5.2.3) 

AAjj - (Ui + UuHVi-V^+dSi) 

NOTE: The formula assumes negligible etror due to changes in the flight track vector from 


iteration (i-1 ) to (i). 


5.3.5 Noise Estimate on Sideline 


Multiple sideline noise estimates are 


included with the noise contour computer program. The 


observer locations for these noise numbers are on 


the sidelines in the (U, V, Z) coordinate system. 


as shown in figure 77. The sideline distances can he specified by the user; the default values are 
1 .0 m. 1 52.4 nf(500 feet), 463.3 m. If any or all of the values r. ed to be changed, they may 


user input. 


Calculation: 


let R; 2 = SD 2 + Z E 2 


a j = arcos(SD/Rj) 
where SD is a set of sideline distances. 

A three-dimensional interpolation on at (EPPj, log Rj, a j) yields the specified 


estimates. 


Boeing Commercial Airplane Company 
P.O. Box 3707 

Seattle, Washington, July 1973 



appendix a 

THEORETICAL GROUND REFLECTION PREDICTION PROCEDURE 


BACKGROUND (PLANE WAVES) 

An acoustic termS fo ' SpheriCa ' 

^ -»*- <*- ^zztzz:z?: ztTz: 

negligible when compared to At even ntultip.es of ISO*. 

:zr "foccuts. This 

or vice-versa. The following is an analysis directed at solvtug this proble . 

Figure A, shows the geometry of the " ^ “p^"^ 

„ is given by the acoustic eguivaien, of 

Snell’s law. 


KqCOS ‘'q-Kj cos ‘'I 


where 


K 0 = 


2fff/C 0 (wave number in air) 


Kj 2*f/C i (wave number in ground) 


In general, K, is complen (non-uniform wave,, but ^ "^i"' 

wave cannot attenuate in the Y-direct,on for both medta (an and grou 

neglected 


* of K affects only the transmitted signal and not the reflected 
real 'part of K, will be considered in wha, follows Spherical 

propagation will be treated later. 



I 


Therefore, 


K, / K q ~ C Q / C-j 


where Cq and C j are speeds of sound in the air and ground respectively. 


Continuing with the solution for sin , yields 


2 ^2 


slnV, = 1 - (W cos 

If Ko/Kj is greater than unity, as is most likely the case for the ground, there is a critical angle, p c , 
in which sin^V| can become negative. That is, 


( Oj FOR v o ^ v c | 
±L FOR ifc < *> c J 


where 


i£ = COS” 1 ( Kj/ K 0 ) 

«,= 1 1 - ( K 0 / ) 2 COS 2 v Q 


The choice of sign v/hen sin v, is imaginary depends upon the convention used to denote phase 
delay in the velocity potentials (0) and impedance (Zj), and in the boundary condition 
transmitted signal (0 j) should vanish as X approaches infinity. 

The convention used here is 

0j = Incident signal 


= * 0 exp[-i K q (X sin v Q + Y cos* 0 )] 


0^ = reflected signal 


0 r exp[-1 K 0 (-X sin v Q + Y cosp q /] 


ansmitted signal 


- 0 fj T exp[=1 K 1 (X sin^ + Y cos# ( )l 

F = plane-wave reflection coefficient 
{Z^/Zq) - (sin / s1n* 0 ) 


yZp) + (sin^ / s 1 n «* 0 j 


y t- V 


T = plane-wave transmission coefficient 

2(VZ 0 ) 

= 1 + r s (Z 1 /Z 0 ) + (sWj 1 sin i^) 


whe „ v 0 is less than - c , we see that* T wilt vanish as X approaches infinity, if and only if; 

= -1 « 1 and cos s ^ * (tyV C0S *o 


where (a,, 0j) are positive real numbers. 


That is: 


0 


T 


lim 0 _ 

x -*• oo T 


0 Q T exp[-1 (“1 ®i x + ^1 
0 q T exp'-Kj a, X) exp(-1 K, Y) 

*H m oo exp(-K^ £«-j X) [0 O T exp(-i K ] fi } Y)] 
0 


These relationships where » T represents a non-U., fonn. ^"^TZTc 
is worth notin, that there is no flow of real acoushc powe » ^ dw ^ . e 

power How in the Y ’ number, ^hen the reflection coefficient T is exp(i 6); 

"reTe^igna. » R is coal in magnitude to the incident s*nai*„ bu, suffers a Phase shtf. 

given by 

9 = 2 tan' 1 !", Z„ / < z l sir,,, o )] 

r;M=s-=: 

critical angle is included in the analysis. 

The composite signal of the direct and reflected sound is then given by the velocity potential 
solution for plane waves as 


0 C ° % + * R = V 1 + 7 V 

- 0 O [1 + r exp(~1 K Q AP)] 


(A3) 
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where P D is .he signal for the direct path end r is the reflection coefficient given in equation (A2>. 

For acoustically hard surfaces like water or concrete, I Z,/Z 0 | may be as large as 3 « I0 3 , *> 
that the reflection coefficient is approximately unity, except for very small angles a u 0 . When v 
approaches aero (grating incidence), the T approaches minus one regardless of the rmpedance 
values. Since the difference in path length, AP, also approaches aero under these c0n< "“ 0 " s - * 

erroneous conclusion given b, equation (A3) is that the observed signal. P c 
contradiction is solved by noting that uniform acoustic plane waves just do not exist. A m re 

, {5 cVvnu/q that the wave fronts are tent in the vicinity of the 
detailed analysis in leferences © and © shows that tne wave irom* <u 

ground for a “dissipative” reflecting plane (see figure *3). 

Since we are interested in the reflected signal, not the transmitted signal, further analysis of 
this 1 layer phenomenon can seive only academic interest. Its presentation was ,0 point om 
™ reason why equation (A3) fails for propagation at grazing incidence. Another approach follows 
whicZisidem acoustic waves with spherical propagation. I. does no, have this scanty. 


SPHERICAL PROPAGATION 

Rudnick (ref. 16) showed that equation (A3) can be replaced by 

*C = # 0 | 1 + r exp('1 K 0 AP)j (A4) 


where 


r ' = effective reflection coefficient for 
spherical propagation 

AP v 

= [r+ (i -r) F(w)] / (i + -j) 

r = plane-wave reflection coefficient 
F(W) = the "boundary loss factor" 

~ 1 + 12 4$ exp(- d) J exp(-Z ) dZ 

= 1 » fF Z exp(Z 2 ) Erfc(Z) with Z 53 -1 VB 
W = [i 2 K P sin 2 **,] / [ 7 ^* (1 -r) COS V Q ] 2 

O * **Q 
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Equation (A4) applies to & point source, homogeneous media (air and ground) and a 

•smooth/infinite/reflecting plane with complex acoustic wave impedance. When the observer ge 
smooth/infimte/reiiecting piane v annroaches zero If the angle of incidence, v 0 , 

sufficiently far away from the sound source, &P/P approaches zero. 11 rnc ^ & o 

approaches 90° under these conditions, then equation (A4) become idem.cal to equation (A3). 

L "pLe wave apptoaimatlon is aa *»d as M for aphencal “ 

irirciaft/observer geometry cannot be restricted to jost this case and equatton <A4) re be« 

Z2Z fo, describing rite ground reflection phenomenon and it has therefore been 

incorporated into the extrapolation methods used in the computer program. 


For the purpose of evaluating the “boundary loss factor ” the following expansions are given. 


When Iwlis less than ten, then 


W K 

( 2 K+ 1 ) K! 


(A5.1) 


When I W| is greater than or equal to ten, then 

F(W) » 9 (W) / [W + 9 i")] 


with 9(H) * -[0.5 + 

W - H .5 “ 4 


(A 5 . 2 ) 


VJ 6 *'5 Tb. 803 

U - 4 .'269.] 


Due to truncation error by the computer (not enough significant digits), equation (A5.1) is not 

-ST- - - ->r— . W L-J r t :=== 


ran rCCiAS 2 , is Stabie for values of |w|grea.er than .0. it was formed by rimreforming 
the Tavior series of Erfc(Z) into Gaussian continued fraction (« f - W and truncating a er ive 

terms. The maximum absolute error of this approximation is (1.1 x 10 ). 


bandwidth effects 


I- - f**" 


So far the anaiysis has considered .mi, a simple harmonic source. Wha, -uriw . en^ueere are 
really interested in is how ground reflection affects Sound Pressure U vel ipcctra «PLS» a 
measured with finite bandwidth equipment. Since detection equipment sums the signals with 


■// : i v 


I 


t 


frequencies contained in the bandwidth of the filters, it is necessary to integrate equation (A4) to 
determine the bandwidth effect. 


ASPL s 10 log 10 



* 


I 


I 


t 


where (•) denotes complex conjugate and (f L , denote the cutoff frequencies for the filters. 

If the values for Zj/Zq and K,/K 0 do not vary erratically over the frequency limits of 
integration, this integral can be approximated by 


r 9 Sln($.,)] 

ASPL * 10 Log 10 ll + A + 2 A cos(S 2 -8) — jj— J 


(A6) 


where 

A *|r| and0=arg(r f ) 

S, -*AP(fu' f L^ c 0 

= n AP(fu + f^/CQ 

and the values for r are evaluated at the geometric-mean-frequencies tor the filters. 


t 
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OTHER EFFECTS 


There are other effects which affect noise measurements in the vicinity of a reflecting ground 
une: wind and temperature gradients: noise source distribution : and scattering of the sound by an 
ousticatly rough ground plane, etc. These other effects have rot been adequately quantified 
is time for incorporation into mathematical terms. 


I 
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GROUP INDEX 727 Retrofit Feasibility IVogram 

SUBJECT Theoretical Ejector Performance Parameter* 

Utilizing Simulated JT8D-9 Engine Condition* 


In order to moke certain ejector design decisions, the 
hell a f rc ,h th«l° f low *J[o^rt7« ©'Xormined by a 

the figures: 


Where: 


A M /A p 

$ 

M E 
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V T T- 

Ps/Pc» 

P !%> 

v 
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Fig. 7 

Fig. »5 

F*g. 23 
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Fig 8 
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These corves are fo provide trend data only and ore not to be util ized 
estimating the performance of o specific ejector-suppressor conf.gurat.on. 

It should be noted that the primary flow Mach Number ot Station © 15 no 
the fully expanded Mach Number. 

Unfortunately, this ejector program assumes con. K !s*« mixing and -Joes not 
provide the ejector flow properties between the entrance a..w •«■» <* the 
shroud. A mixing program is being formulated by the fcopulsion Researcn 
group to provide flow properties as a function of axiol displacement. T e 
mixing program will probably be ready for checkout in December 1971. 


Prepared by 



Approve 


; Vjp -- 

>v ©a DV 


E. Tjonneland 


Attachment: Table I, Figures 1-32 
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TABLE I 



EJECTOR LENGTH-DIAMETER RATIO, L/D = 2.5 
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